Fluctuating Stripes in Strongly Correlated Electron Systems and the 
Nematic-Smectic Quantum Phase Transition 
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We discuss the quantum phase transition between a quantum nematic metallic state to an electron 
metallic smectic state in terms of an order-parameter theory coupled to fermionic quasiparticles. 
Both commensurate and incommensurate smectic (or stripe) cases are studied. Close to the quantum 
critical point (QCP), the spectrum of fluctuations of the nematic phase has low-energy "fluctuating 
stripes". We study the quantum critical behavior and find evidence that, contrary to the classical 
case, the gauge-type of coupling between the nematic and smectic is irrelevant at this QCP. The 
collective modes of the electron smectic (or stripe) phase are also investigated. The effects of the 
low-energy bosonic modes on the fermionic quasiparticles are studied perturbatively, for both a 
model with full rotational symmetry and for a system with an underlying lattice, which has a 
discrete point group symmetry. We find that at the nematic-smectic critical point, due to the 
critical smectic fluctuations, the dynamics of the fermionic quasiparticles near several points on the 
Fermi surface, around which it is reconstructed, are not governed by a Landau Fermi liquid theory. 
On the other hand, the quasiparticles in the smectic phase exhibit Fermi liquid behavior. We also 
present a detailed analysis of the dynamical susceptibilities in the electron nematic phase close to 
this QCP (the fluctuating stripe regime) and in the electronic smectic phase. 

PACS numbers: 71.10.Hf, 71.45.Lr, 71. 10. Ay 



I. INTRODUCTION 

The discovery of the high temperature superconduc- 
tors in the quasi-two-dimensional copper-oxide materials 
in the late 1980s, and of novel correlated phases in other 
complex oxides, has brought to the forefront the prob- 
lem of the physics of strongly correlated electron systems. 
To this date the understanding of the behavior of these 
systems remains one of the main open and challenging 
problems in condensed matter physics. The central co- 
nundrum in this field is the fact that these strongly cou- 
pled electron systems are best regarded as doped Mott 
insulators for which both the band theory of metals and 
the Landau theory of the Fermi liquid (FL) fail. 

One characteristic feature of the physics of doped 
Mott insulators is their inherent tendency to electronic 
phase separation, frustrated by the effects of Coulomb 
interactions, [l], 0] The ground states resulting from these 
competing tendencies typically break the translation in- 
variance and/or the point group symmetry of the un- 
derlying lattice. From a symmetry point of view, the 
ground states of doped Mott insulators are charge- 
ordered phases, which share many similarities with clas- 
sical liquid crystals, and should be regarded as electronic 
liquid crystal phases. Q However, unlike classical liquid 
crystals, electronic liquid crystals are strongly quantum 
mechanical states whose transport properties range from 
insulating to metallic and even superconducting. In con- 
trast with classical liquid crystals, whose ordered phases 
represent the spontaneous breaking of the continuous 
translation and rotational symmetry of space 0, Q, the 
electronic liquid crystal phases of strongly correlated sys- 
tems are sensitive to the effects of the underlying lattice 



and the symmetry breaking patterns involve the point 
and space groups, as well as to disorder. More complex 
ordered states, involving simultaneously charge and spin 
degrees of freedom, may also arise. [1] 

The sequence of quantum phase transitions described 
above, electron crystal — ► smectic (stripe) — ► nematic — > 
isotropic fluid, representing the progressive restoration 
of symmetry, is natural from a strong correlation per- 
spective. Indeed, the electron crystal state(s) are natu- 
rally insulating (much as in the case of a Wigner crys- 
tal), the smectic or stripe phases are either anisotropic 
metals or superconductors, and the charged isotropic 
fluids are either metallic or superconducting. While 
the isotropic metallic phase is essentially a FL (albeit 
with strongly renormalized parameters), the nematic and 
smectic metallic phases have a strong tendency to show 
non-FL character. Indeed, much of the theoretical de- 
scription of the stripe or smectic phases is usually based 
on a quasi-one-dimensional analysis, which makes ex- 
plicit use of this strong correlation physics. Such ap- 
proaches give a good description of this state deep in- 
side this phase and at energies high compared to a 
"dimensional crossover" scale below which the state is 
fully two-dimensional (and strongly anisotropic) 0, H, 
S 011 EH, E2, EH- Stripe phases (insulating, metal- 
lic, and superconducting) have been found in mean-field 
studies of generalized two-dimensional Hubbard and t-J 

models [3 M El, Ei EE m M M M M H HE m . 

The same pattern of quantum phase transitions can 
also be considered in reverse order, with a weak coupling 
perspective, as a sequence of symmetry breaking phase 
transitions beginning from the isotropic metal: FL — > 
electron nematic — > electron smectic — > insulating elec- 



2 



tron crystal. In this case, one begins with a uniform 
isotropic metal, well described at low energies by the Lan- 
dau theory of the FL, with well-defined quasiparticles and 
a Fermi surface (FS), and considers possible instabilities 
of the isotropic fluid into a nematic (or hexatic and other 
such states), as well as phase transitions into various pos- 
sible charge-density-wave (CDW) phases. The unidirec- 
tional CDW-ordered states are the weak coupling analog 
of the smectic (or stripe) phases, and have the same or- 
der parameters as they break the same symmetries. The 
main difference between a CDW and a smectic resides in 
the fact that while the CDW arises as a weak coupling 
(infinitesimal) instability of a FL in which parts of the FS 
are gapped[27| (which requires the existence of a FS with 
sharp quasiparticles), the stripe phases do not require 
such description. While a CDW phase at high energies is 
essentially a FL, the high-energy regime of a stripe phase 
is a quasi-one-dimensional Luttinger liquid. [8|, [9| A direct 
quantum phase transition from a FL to a CDW phase is, 
naturally, possible and this quantum phase transition has 
been studied in some detail, [28, 29] as well as to a metal- 



lic spin-density wave (SDW) [3d. l31| 



The weak coupling description of an electron nematic 
phase uses a Pomeranchuk instability of a Fermi liq- 
uid state [32|. Oganesyan, Kivelson, and Fradkin[33j 
showed that the nematic quantum phase transition is 
a quadrupolar instability of the FS, and gave a char- 
acterization of the properties of the nematic Fermi fluid 
in a continuum model. An electron nematic quantum 
phase transition has also been found in lattice models [34], 
l35l [3^ . which show, however, a strong tendenc y to ex- 
hibit a first-order quantum phase transition [37l l38l. [39l|. 
Pomeranchuk instabilities in the Landau theory of the 
FL have also shown the existence of an electron ne- 
matic transition [40l. |4lj |. Perturbative renormalization 
group analysis of the stability of the FL in Hubbard-type 
models 42], as well as high-temperature expansions [43], 
has also shown that in such models there is a strong ten- 
dency to a nematic state. An electron nematic state was 
shown to be the exact ground state in the strong cou- 
pling limit of the Emery model of the copper oxides at 
low hole doping [3]. 

The upshot of the work on the electron nematic quan- 
tum phase transition is that, at the QCP (if the transi- 
tion is continuous) and in the nematic phase (in the con- 
tinuum) the electron quasiparticle essentially no longer 
exists as an asymptotically stable state at low energies, 
except along symmetry determined directions in the or- 
dered phase. A full solution of this QCP by bosonization 
methods has confirmed these results, which were gleaned 
from mean-field theory, and have also provided stron g ey - 
idence for local quantum criticality at this QCP [45ll46|]. 

In this paper we will be interested in the quantum 
phase transition from an electron nematic phase to a 
charge stripe phase, a unidirectional CDW. For simplic- 
ity we will not consider here the spin channel, which plays 
an important role in many systems. We will only con- 
sider the simpler case of unidirectional order. Extensions 



to the more general case of multidirectional order are 
straightforward. Here we develop a quantum-mechanical 
version of the nematic-smectic transition in a metallic 
system. This is a quantum-mechanical version of the 
McMillan-deGennes theory for the quantum phase tran- 
sition from a metallic nematic phase to a metallic smectic 
(or CDW) phase. The construction of such a generaliza- 
tion of the McMillan-deGennes theory is the main pur- 
pose of this paper. 

As it is discussed in detail in subsequent sections, here 
we will follow the "weak-coupling" sequence of quantum 
phase transitions described above, beginning with the 
transition from a FL to an electron nematic, and from 
the latter to a stripe or unidirectional CDW state. The 
main advantages of this approach are that it allows to 
address the fate of the electronic quasiparticles and non- 
Fermi liquid behaviors as the correlations that give rise 
to these electronic liquid crystal phases develop, as well 
as to study the quantum critical behavior following the 
standard Hertz-Millis approach [i?], O [Zi|. However, 
the main disadvantage is that this approach does not do 
justice to the physics of strong correlation. For this rea- 
son, in spite of the important insights that are gained 
through this line of analysis, this approach cannot ex- 
plain the physics of the "strange metal" regime observed 
in the "normal state" of high T c superconductors where 
non-Fermi liquid effects are widely reported. To do that 
would require studying this problem as a sequence of 
quantum melting transitions. An important first step in 
this direction has been made by Cvetkovic and cowork- 
ers [lii HH, H3] who have studied a purely bosonic model 
of such quantum melting. The inclusion of fermionic de- 
grees of freedom in this strong coupling approach is an 
interesting but challenging open problem. 

We have both conceptual and phenomenological mo- 
tivations for considering this problem. At the concep- 
tual level the main question is to develop a theory of 
the quantum critical behavior at the electron nematic- 
smectic phase transition, and of the low-energy physics 
of both phases near quantum criticality. Although the 
static properties are the same as in the classical the- 
ory (as required by symmetry) the quantum dynamics 
changes the physics substantially. Thus, physical prop- 
erties, which determine the transport properties and the 
fermion spectral function, cannot be gleaned from the 
classical problem. Provided that the quantum phase 
transition is continuous or, at most weakly first order, 
the low-energy fluctuations in one phase (say the ne- 
matic metal) must reflect the character of the nearby or- 
dered stripe phase. In other words, under these assump- 
tions, as the quantum phase transition is approached the 
metallic nematic phase behaves as a state with "fluctu- 
ating stripes" . The ample experimental evidence in high 
temperature superconductors for "fluctuating stripe or- 
der" should be interpreted instead as evidence of a ne- 
matic phase proximate to a quantum phase transition to 
a stripe (or smectic)-ordered state [53| . 
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II. SUMMARY OF RESULTS 

In this work we follow a phenomenological approach 
to study the quantum phase transition between an elec- 
tronic nematic state and electronic smectic state. We 
postulate the existence of both an electron nematic and 
a smectic phases with a possible direct phase transition 
between them. This physics will be represented by an 
effective field theory involving the nematic and CDW or- 
der parameters. The static part of the effective action 
of the order-parameter theory has the same form as in 
the classical theory of the nematic-smectic transition, the 
McMillan-deGennes theory. We will assume that aside 
from the effects of the coupling to the fermionic quasi- 
particles, this effective field theory is analytic in the or- 
der parameters and their derivatives as this dependence 
is determined by local physics. As shown below, this as- 
sumption implies a dynamical quantum critical exponent 

2=1. 



The fermionic quasiparticles couple to the nematic 
and smectic (CDW) order parameters in their natural 
symmetry-dictated way. The fermions will be assumed 
to be a normal FL, with well-defined quasiparticles and 
a FS. Thus, we will not attempt to explain why the phase 
transition exists, which requires a microscopic theory, but 
rather describe its character. One of our most important 
results is that this theory gives a description of a phase 
with fluctuating stripe (smectic) order, of much interest 
in current experiments. The effective theory that we con- 
sider also allows for a possible direct transition between 
the normal and isotropic FL state and a CDW phase, 
without going through an intermediate nematic phase, 
as in the direct transition between a FL and a CDW 
state, discussed by Altshuler, Ioffe, and Millis [28| . Thus, 
the theory we present here actually describes the behav- 
ior of a FL in the vicinity of a possible bicritical point 
which, as we shall see, is not directly accessible. 
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TABLE I: Summary of results. See the text for a detailed explanation. 



The main results of our theory are summarized in Ta- 
ble HI In Sec. IIIII we discuss the current experimental 
status of electronic liquid crystal phases in a number of 
different materials. In Sec. IIVI we set up the order pa- 
rameter theory for the electronic liquid crystal phases 
based on symmetry and analyticity. The static part of 
this phenomenological theory is (as it should be) similar 
to its classical counterpart, but we add proper dynamics 
to describe the quantum fluctuations. We next couple 
the order parameter theory to the fermionic quasipar- 
ticles, in Sec. |Vj The coupling between the fermionic 
quasiparticles and the order parameters is completely de- 
termined by symmetry. This is a standard approach to 
study quantum phase transitions in metallic systems [4^ | . 
It is a consistent scheme for the study of the quantum 
phase transition provided the effective dimension d + z is 
close to 4 (here d is the dimensionality of space). Several 
different non-analytic dependences on the order param- 
eters in the effective action appear as a consequence of 
their coupling to the fermions. We show that these non- 



analytic dynamical terms dominate over the dynamics 
prescribed phenomenologically. Hence, the dynamics of 
fermionic liquid crystal phase is very different from that 
of the simple phenomenological theory. We present a 
detailed analysis of the behavior of the dynamical sus- 
ceptibilities in both phases and at the QCP. 

The nematic-smectic QCP is studied in Sec. I VII In 
classical liquid crystals, the Goldstone mode of the ne- 
matic phase plays a very important role at the nematic- 
smectic transition. There, this relevant coupling drives 
the transition weakly first order through a fluctuation- 
induced first order transition 54]. However, in the case 
of the electronic liquid crystals, we find that the cou- 
pling between the nematic Goldstone mode and the smec- 
tic field is actually irrelevant at the electronic nematic- 
smectic QCP. Therefore, these two modes can be treated 
separately, as they are weakly coupled to each other. Sev- 
eral different nematic-smectic critical theories are stud- 
ied, depending on the relation between the magnitude 
of the ordering wave vector of the CDW, Qs, and the 
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Fermi wave vector, kp. For Qs < 2kp (Fig. HJa)), we 
find that the critical smectic field has a dynamic critical 
exponent z — 2, which will result in a C ~ T contribution 
to the heat capacity. This is a correction to the conven- 
tional linear T behavior of Fermi liquids. These quan- 
tum fluctuations lead to the existence of four points on 
the FS where the assumptions of FL theory are violated 
(Fig. UK a)). At these points the imaginary part of the 
fermion self-energy correction T,"(kp,uj) ~ M 1 / 2 . For 
Qs = 2kp (Fig. EJb)), the system exhibits anisotropic 
scaling: [q x ] = 1, [q y ] = 2 and [u>] = 3 for the incom- 
mensurate CDW, while [q x ] — 1, [q y ] — 2 and [u>] = 2 
for the commensurate case. Besides, a non-analytic $ 5 / 2 
term, where $ is the smectic order parameter, is gen- 
erated in the action of the low-energy effective theory. 
This non-analytic term is relevant under the renormaliza- 
tion group (RG) for the incommensurate case, suggesting 
a weak, fluctuation-induced, first-order transition. This 
coupling is irrelevant in the commensurate case. Here we 
also find two points on the FS (Fig. |U[b)), where the 
system has marginal FL behavior, with a quasiparticle 
scattering rate E"(fcp,w) ~ \lo\, and a low temperature 
correction to the heat capacity C ~ T 3 / 2 , which is sub- 
leading. We also consider the special case of a CDW 
caused by a nearly nested FS, for which we find that 
the low-temperature heat capacity correction C ~ T 4 / 3 , 
which is also subleading, and the fermions form a FL, 
with Y,"{k F ,ui) ~ H 13 / 12 . We also calculated the dy- 
namic CDW susceptibility Xs{Q.i u ) f° r both cases. The 
Qs > 2kp case will not be discussed here. In the pres- 
ence of a lattice this case is quite trivial (see Sec. IVIj) 
while for it to occur in a continuum system, where it is 
non-trivial, requires unphysical assumptions. 

The smectic phase is discussed in Sec. IVIII In the 
smectic phase the anisotropic scaling associated with the 
Goldstone fluctuations are [q x ] = 1, [q y ] = [u>] = 2. We 
find that the low-temperature heat capacity correction 
C ~ T 3 / 2 , which is also subleading. The quasiparticle 
scattering rate in this case is T,"(kp,u>) ~ w 2 log|o;| for 
much of the FS while E"(/cf, w) ~ |w| 3 / 2 at the two spe- 
cial points where the Fermi velocity is parallel to the 
ordering wave vector. Thus, in this case fermions behave 
as a FL. We also calculated both the longitudinal and 
transverse dynamic CDW susceptibilities in the smectic 
phase. 

Lattice effects are also discussed. For the case of an 
incommensurate smectic phase, we show that there is 
an unpinned smectic phase close to the nematic-smectic 
critical point. In this phase, the smectic Goldstone mode 
has a dynamic critical exponent z = 1 and the system 
is a FL with E"(w) ~ w 2 log|w| at most of the FS and 
E"(o;) ~ u! 2 at some special point on the FS described 
below. Due to the unpinned smectic ordering, the system 
receives a correction to the low-temperature heat capac- 
ity C ~ T 2 , and we also computed the dynamic trans- 
verse CDW susceptibility. Deep into the smectic phase, 
an incommensurate CDW may be pinned down by lattice 
distortion. As expected, the fermions in a pinned smectic 



are in a conventional FL state. 

In Sec. I Villi we present a brief discussion of the role of 
thermal fluctuations for these phases and of the classical- 
to-quantum crossovers. We conclude with a summary of 
our main results and a discussion of open questions in 
Sec. IIX1 Details of the calculations are presented in 
several appendices. In Appendix [A] we discuss the ten- 
sor structure of the order parameters. In Appendix [B] 
we present details of the nematic-smectic QCP for the 
case Qs < 2kp, while the nonanalytic terms induced for 
the Qs = 2kp case are presented in Appendix [U] In 
Appendix [D] we present details of the calculation of the 
spectrum of Goldstone modes in the smectic phase. In 
Appendix [E] we summarize the random phase approxi- 
mation (RPA) calculation of the fermion self-energy at 
the nematic-smectic QCP and in the smectic phase. 



III. EXPERIMENTAL STATUS OF 
ELECTRONIC LIQUID CRYSTAL PHASES 

During the past decade or so experimental evidence has 
been mounting of the existence of electronic liquid crystal 
phases in a variety of strongly correlated (as well as not as 
strongly correlated) electronic systems. We will be par- 
ticularly interested in the experiments in the copper ox- 
ide high temperature superconductors, in the ruthenate 
materials (notably S^R^Oy), and in two-dimensional 
electron gases (2DEG) in large magnetic fields. However, 
as we will discuss below, our results are also relevant to 
more conventional CDW systems such as the quasi-two- 
dimensional dichalcogenides. 



A. High temperature superconductors 

In addition to high temperature superconductivity, the 
copper oxide materials display a strong tendency to have 
charge-ordered states, such as stripes. The relation be- 
tween charge ordered states (55|, as well as other proposed 
ordered states[56|,[57|], and the mechanism(s) of high tem- 
perature superconductivity is a subject of intense current 
research. It is not, however, the focus of this paper. 

Stripe phases have been extensively investigated in 
high temperature superconductors and detailed and re- 
cent reviews are available on this subject [H, HI]. Stripe 
phases in high temperature superconductors have unidi- 
rectional order in both spin and charge (although not 
always) and it is typically incommensurate. In gen- 
eral the detected stripe order (by low-energy inelastic 
neutron scattering) in La2_, T Sr 7 ,Cu04, La9 -. r Ba, T CuQ4 
and YBa2Cu306+ y (see Refs.[53[ and [58[ and refer- 
ences therein) is not static but "fluctuating". As em- 
phasized in Ref . [53l] . "fluctuating order" means that 
there is no true long range unidirectional order. In- 
stead, the system is in a (quantum) disordered phase, 
very close to a quantum phase transition to such an or- 
dered phase, with very low-energy fluctuations that re- 
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veal the character of the proximate ordered state. On 
the other hand, in L^-zBa^CuC^ near x = 1/8 (and 
in Lai.6-a;Ndo.4Sr a ;Cu04 also near x = 1/8), the order 
detected by elastic neutron scattering [59(, and resonant 
x-ray scattering in La2_ a; Ba 2 ;Cu04 60] also near x = 1/8, 
becomes true long-range static order. 

In the case of La2- £C Sr a: Cu04, away from x = 1/8, and 
particularly on the more underdoped side, the in-plane 
resistivity has a considerable temperature-dependent 
anisotropy (6ll]. which has been interpreted as an indi- 
cation of electronic nematic order. From these experi- 
ments it has been suggested that this phase be identified 
as an electron nematic [6l|. The same series of experi- 
ments also showed that very underdoped YBa 2 Cu306+j / 
is an electron nematic as well. The most striking evidence 
for electronic nematic order in high temperature super- 
conductors are the recent neutron scattering experiments 
in YBa2Cu306+j, at y = 6.45[62|]. In particular, the 
temperature-dependent anisotropy of the inelastic neu- 
tron scattering in YBa2Cu3 0g+y shows that there is a 
critical temperature for nematic order (with T c ~ 150K) 
where the inelastic neutron peaks also become incom- 
mensurate. Similar effects were reported by the same 
group [63| at higher doping levels (y ~ 6.6) who ob- 
served that the nematic signal was decreasing in strength 
suggesting the existence of a nematic-isotropic quantum 
phase transition closer to optimal doping. Fluctuating 
stripe order in underdoped YBa 2 Cu306+ y has been de- 
tected earlier on in inelastic neutron scattering experi- 
ments [&3.l65j| which, in hindsight, can be reinterpreted as 
evidence for nematic order. However, as doping increases 
the strength of the temperature-independent anisotropic 
background, due to the increased orthorhombicity of the 
crystal, also increases thus making this phase transition 
difficult to observe. 

Recent inelastic neutron scattering experiments have 
found similar effects in La 2 _ a; Sr a ;Cu04 materials where 
fluctuating stripes where in fact first discovered [66] . Mat- 
suda et al [67| have given qualitatively similar evidence 
for nematic order in underdoped La^zSr^CuC^ (x = 
0.05) which was known to have "fluctuating diagonal 
stripes" . In the same doping range it has also been 
found by resonant x-ray scattering experiments that 5% 
Zn doping stabilizes a static diagonal stripe-ordered state 
with a very long persistence length which sets in at quite 
high temperatures (68|. 

These recent results strongly suggest that the experi- 
ments that had previously identified the high tempera- 
ture superconductors as having "fluctuating stripe order" 
(both inside and outside the superconducting phase) were 
most likely detecting an electronic nematic phase, quite 
close to a state with long-range stripe (smectic) order. 
In all cases the background anisotropy (due to the or- 
thorhombic distortion of the crystal structure) acts as a 
symmetry breaking field that couples linearly to the ne- 
matic order, thus rounding the putative thermodynamic 
transition to a state with spontaneously broken point 
group symmetry. These effects are much more appar- 



ent at low doping where the crystal orthorhombicity is 
significantly weaker. 

The nature of the fluctuating spin order changes sub- 
stantially as a function of doping: in the very underdoped 
systems there is no spin gap while inside much of the su- 
perconducting dome there is a finite spin gap. In fact in 
La 2 _ 2: Ba 2 ;Cu04 at x — 1/8 there is strong evidence for 
a complex stripe-ordered state which combines charge, 
spin and superconducting order (69l Fro] ]. These experi- 
ments have also established that static long-range stripe 
charge and spin orders do not have the same critical tem- 
perature, with static charge order having a higher T c . 

An important caveat to our analysis is that in doped 
systems there is always quenched disorder, and has dif- 
ferent degrees of short range "organization" in differ- 
ent high temperature superconductors. Since disorder 
also couples linearly to the charge order parameters 
it ultimately also rounds the transitions and renders 
the system to a glassy state (as noted in Refs.fl, HU). 
Such effects are evident in scanning tunneling microscopy 
(STM) experiments in Bi 2 Sr 2 CaCu208+<5 which revealed 
that the high-energy (local) behavior of the high tem- 
perature superconductors has charge order and it is 
glassy 0,0, S Izllzj. 

Finally, we note that in the recently discovered iron 
pnictides based family of high temperature supercon- 
ductors, such as La (Oi-^F^FeAs [75J, |76j, a unidirec- 
tional spin-density wave has been found. It has been 
suggested[77| that the undoped system LaOFeAs may 
have a high-temperature nematic phase and that quan- 
tum phase transitions also occur as a function of fluorine 
doping[78]. This suggests that many of the ideas and re- 
sults that we present here may be relevant to these still 
poorly understood materials. 



B. Other complex oxides 

The existence of stripe-ordered phases is well estab- 
lished in other complex oxide materials, particularly the 
manganites and the nickelates. In general, these ma- 
terials tend to be "less quantum mechanical" than the 
cuprates in that they are typically insulating (although 
with interesting magnetic properties) and the observed 
charge-ordered phases are very robust. These materi- 
als typically have larger electron-phonon interactions and 
electronic correlations are comparatively less dominant 
in their physics. For this reason they tend to be "more 
classical" and less prone to quantum phase transitions. 
However, at least at the classical level, many of the issues 
we discussed above, such as the role of phase separation 
and Coulomb interactions, also play a key role 79]. The 
thermal melting of a stripe state to a nematic has been 
seen in the manganite material Bi^Ca^MnOsfSO]. 
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C. Ruthenates 

Recent magneto-transport experiments in the quasi- 
two-dimensional bilayer ruthenate Sr3Ru207 by the St. 
Andrews group [81] have given strong evidence of a strong 
temperature-dependent in-plane transport anisotropy in 
these materials at low temperatures T < 800 mK and 
for a window of perpendicular magnetic fields around 
7.5 Tesla. These experiments provide strong evidence 
that the system is in an electronic ncmatic phase in 
that range of magnetic fields [SjJ, ]82j . The electronic ne- 
matic phase appears to have preempted a metamagnetic 
QCP in the same range of magnetic fields |j| 184L I85L [8?| . 
This suggests that proximity to phase separation may 
be a possible microscopic mechanism to trigger such 
quantum phase transitions, consistent with recent ideas 
on the role of Coulomb-frustrated phase separation in 
2EDGs[13,[8l. 



D. 2DEGs in large magnetic fields 

To this date, the best documented electron nematic 
state is the anisotropic compressible state observed in 
2DEGs in large magnetic fields near the middle of a Lan- 
dau level, with Landau index N > 2[H H EEO, H2]. 
In ultrahigh-mobility samples of a 2DEG in AlAs-GaAs 
heterostructures, transport experiments in the second 
Landau level (and above) near the center of the Lan- 
dau level show a pronounced anisotropy of the longi- 
tudinal resistance rising sharply below T ~ 80 mK, 
with an anisotropy that increases by orders of magni- 
tude as the temperature is lowered. These experiments 
were originally interpreted as evidence for a quantum 
Hall smectic (stripe)_phase [H, [§! S, Fur- 
ther experiments jiM liM llOOj did not show any evidence 
of pinning of this putative unidirectional CDW as the 
I —V curves are strictly linear at low bias and no broad- 
band noise was detected. In contrast, extremely sharp 
threshold electric fields and broadband noise in transport 
was observed in a nearby reentrant integer quantum Hall 
phase, suggesting a crystallized electronic state. These 
facts, together with a detailed analysis of the experimen- 
tal data, suggested that t he compressible stat e is in an 
electron nematic phased Eml ll02L 1 103L Il04l | , which is 
better understood as a quantum melted stripe phase. 



E. Conventional CDW materials 

CDWs have been extensively studied since the mid- 
seventies and ther e are extensive reviews on their 
properties [l05|, Il06j |. From the symmetry point of view 
there is no difference between a CDW and a stripe (or 
electron smectic). The CDW states are usually observed 
in systems which are not particularly strongly corre- 
lated, such as the quasi-one-dimensional and quasi-two- 
dimensional dichalcogenides, and the more recently stud- 



ied tritellurides. These CDW states are reasonably well 
described as FLs which undergo a CDW transition, com- 
mensurate or incomm ensur ate, triggered by a nesting 
condition of the FS [27l . Il07l |. As a result, a part or all of 
the FS is gapped in which case the CDW may or may 
not retain metallic properties. Instead, in a strongly 
correlated stripe state, which has the same symmetry 
breaking pattern, at high energy has Luttinger liquid 
behavior [li, [3. 

What will interest us here is that conventional quasi- 
2D dichalcogenides, the also quasi-2D tritellurides and 
other similar CDW systems can quantum melt as a 
function of pressure in Ti Se9 1 10811 . or by chemic al in - 
tercalation as in Cu^TiSe^lod IllO] and Nb a TaS 2 [m|]. 
Thus, CDW phases in chalcogenides can serve as a weak- 
coupling version of the problem of quantum melting of 
a quantum smectic. Interestingl y, the re is strong experi- 
mental evidence that both TiSe 2 [l0i] and NbzTaS 2 |lIl| 
do not melt directly to an isotropic Fermi fluid but go 
instead through an intermediate phase, possibly hexatic. 
(Cu a: TiSe2 is known to become superconducting [ lQOj ] .) 
Whether or not the intermediate phases are anisotropic 
is not known as no transport data is available in the rel- 
evant regime. 

The case of the CDWs in tritellurides is more di- 
rectly relevant to the theory we present in this pa- 
per. Tritellurides are quasi-2D materials which for a 
broad range of temperatures exhibit a unidirectional 
CDW (i.e. an electronic smectic phase) and whose 
anisot ropi c behavior appears to be primarily of electronic 
origindll [ill El [US Ell- However, the quantum 
melting of this phase has not been observed yet. Theo- 
retical studies have also suggested that it may be possible 
to have a quantum phase transiti on to a state with more 
than one CDW in these materials 1171 ] . 



IV. ORDER-PARAMETER THEORY 

In this section we will construct, using phenomenolog- 
ical arguments, an effective order parameter theory that 
will describe both the electron nematic and the electron 
smectic (or unidirectional CDW) phases. Although by 
symmetry the order parameter theory must be very sim- 
ilar to the ones used in classical liquid crystal phases, 
we will go through the construction of the phenomeno- 
logical theory in some detail for several reasons. In 2D 
the rotation group SO(2) is Abelian which allows for a 
significant simplification of the formulas by using a com- 
plex order parameter for the nematic phase, instead of 
a tensor expressions commonly used for 3D classical liq- 
uid crystals. Proper dynamical terms now need to be 
included to describe the quantum fluctuations at zero- 
temperature. Besides, in order to provide a clear relation 
between this paper and earlier studies of the CDW state 
of fermions, we would like to discuss also the relation 
between the smectic phase and the CDW state. 
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A. The normal-electronic nematic transition 

The nematic order parameter in 2D is a I = 2 repre- 
sentation of the SO(2) rotational group [331 ]. It is defined 
as a symmetric traceless tensor of rank two. 

N= ( Ul1 ni2 V (4.1) 

The 2D rotational group ,90(2) is isomorphic to U(l). 
Hence, we define instead the complex order-parameter 
field N(r, t) 

N(F,t)=n ll (f,t)+in 12 { i r,t), (4.2) 

where r and t are the space and time coordinates. We 
will use this complex order parameter field in this paper 
to take the advantage of the Abelian nature of SO(2). 

The conjugate field is N^(r,t) = nn(f,t) — ini 2 (r, t). 
Under a global rotation by an angle 0, the fields 
N(r, t) and N^r, t) transform, respectively, as N(r, t) — > 
e 2W N(r,t) and N^(r,t) e~ 2ie N^ (r, t) . Hence, N 
and N^ carry the angular momentum quantum numbers 
l z = 2 and l z = —2, respectively. 

This complex order parameter can be generalized eas- 
ily to other angular momentum channels I 2, but not 
to higher dimensions d > 2, since it relies heavily on the 
special property of the 2D rotational group SO(2). In 
higher dimensions, the rotational group will no longer be 
Abelian, so one will need to use the tensor formula as in 
the classical liquid crystal theories. In Appendix [A] for- 
mulas using the complex order parameter are translated 
into the conventional tensor form for comparison. 

The order-parameter field we just defined is invariant 
under spatial-inversion and time-reversal 

PN{r, i)p- 1 = N(-r, t), (4.3) 
TN(f, t)^ 1 = N(f, -t). (4.4) 

In even spatial dimensions, including 2D in which our 
system lives, a chiral transformation is different from a 
space inversion. To change the chirality in 2D, we can re- 
verse the y direction and keep the x direction unchanged. 
Under this chiral transformation, the nematic field will 
be changed into the conjugate field 

CNix^y^C- 1 =N\x,-y,t). (4.5) 

Here, C is the chiral transformation operator. 

The effective action must preserve the symmetries of 
the system, both continuous, as the translational and ro- 
tational symmetries, and discrete, as the time reversal, 
spatial inversion and chiral symmetries. With the as- 
sumption of analyticity, the action must be 

S N = J ArAt (\d t N\ 2 - \V N\ 2 - & N \N\ 2 - u N \N\^ . 

(4.6) 



Here the dynamical term is quadratic in time deriva- 
tives. This is because the term linear in time derivatives 
—iNWtN + h.c. is not allowed by the chiral symmetry. 
It is the imaginary part of N^dtN. It corresponds to a 
pseudoscalar, and is not allowed. 

In 2D, cubic terms in the nematic field N are not al- 
lowed. Hence, if mat > 0, the normal-nematic transition 
is second order, instead of a first-order transition as in the 
3D case d,Q • For un > and A at > 0, the rotational in- 
variant ground state will be stable. When Ajv becomes 
negative, N will develop an expectation value N with 
module yj— An / (2un), which breaks the 5*0(2) rotation 
symmetry. The residual rotational symmetry would be 
Z 2 . The argument of N determines the direction of the 
nematic order parameter. 

The action of Eq. (|4.6p has an internal U(l) symmetry 
associated with the phase of the complex field N, which 
is not physical. By symmetry, terms of the form 

-k (N\d x + id y )N\d x + id v )N + h.c.) (4.7) 

are allowed [33|, Ill8l | . This kind of terms are irrelevant 
at the QCP and in the isotropic phase, which leads to 
the existence of an "emergent" internal U(l) symmetry 
at quantum criticality. But it will be important in the 
nematic phase, as it makes the two Frank constants to 
attain different values. (This effect is formally analogous 
to the role of spin-orbit interactions in the Schrodinger 
equation: in their absence spin is an internal degree of 
freedom.) This emergent symmetry of the normal phase 
and at the critical point is very important for the classical 
normal-nematic transi tion, especially in the study about 
the fluctuation effects [TTg . [l20l Imf . 

B. The electronic nematic phase 

In the nematic phase, the 50(2) rotational symme- 
try is broken. Hence, we expect the fluctuations of the 
amplitude of the nematic order parameter, SN, to corre- 
spond to a massive mode with an energy gap of — 2Aat 
(Aat < 0), and the fluctuations of the phase, (pN, consti- 
tute the gapless Goldstone mode. Without loss of gen- 
erality, throughout this paper, we assume that N is real 
and positive. This state corresponds to a nematic order 
in the main axis direction. In this state, the action of </>at 
is 

S<p N =N 2 J dfdt {(d t (f>N) 2 -Ki(d x( f> N ) 2 - K 2 {d y M 2 ) . 

(4.8) 

where K± = 1 + 2kN and K 2 = 1 — 2kN are the two 
Frank constants. This action is only valid for small ne- 
matic fluctuations. It cannot be used to study topological 
defects of the nematic phase, known as disclinations[l]- 
The field 4>n has dynamic critical exponent z = 1. This 
makes the effective dimension of this system 3, which is 
above the lower critical dimension of the theory d = 2, 
and nematic order will not be destroyed by fluctuations. 
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C. CDW multi-critical point 

The smectic order is a unidirectional CDW, described 
by a single complex order parameter field. If we assume 
analyticity, the effective low-energy theory of the bosonic 
held p can be determined as: 

5 C dw = J drdt(dtp) 2 + S 2 + S 3 + S 4 , (4.9) 

where S2 is the term in the quadratic order of p, S3 and 
S4 are the cubic and quartic terms respectively. 
The term S2 in the momentum space is, 

S2 = ~ J ^f(\k\)p(k, u)p(-k, -w). (4.10) 

The function /(|fc|) has the physical meaning of the in- 
verse of the CDW susceptibility. If we assume the order- 
ing wave- vector of the CDW is Qs (with \Qs\ = Qs its 
magnitude), /(|fc|) will have the form 

/(|fc|) = A CDW + C(\k\ 2 - Q%f + ... (4.11) 

wher e Acdw is the energy gap of the CDW excitations 
[l22j | . and C is a positive constant. 

When the energy gap Acdw decreases to zero, all the 
density wave modes with \k\ = Qs will become soft and 
critical when Acdw = 0. This is very different from an 
ordinary </> 3 - or 4 -theory, where we only need to consider 
one mode (or two modes for a complex field) at small mo- 
mentum. Here, we need to consider all the modes with 
the wave vector k whose magnitude is close to Qs- In 
other words, the point Acdw = is not a critical point 
but a multi-critical point with an infinite number of crit- 
ical modes. Even if a lattice background is present, the 
Acdw = point may still be a multi-critical point of n 
critical modes, if the lattice has a n-fold rotational sym- 
metry for n > 2. For a multi-critical point, higher-order 
terms become important. Without a detailed knowledge 
of these higher-order terms, it is not possible to deter- 
mine whether the transition is first or second order, or 
how many CDWs will be formed in the ordered phase. 

Brazovskii jl22j | studied the classical version of this 
problem, considering only the isotropic interactions. 
Chubukov and co-workers [29j studied the quantum prob- 
lem in a fermionic system in the high-density regime 
where the cubic and quartic terms of p can be ignored. 

In general, depending on the non-Gaussian term s, th e 
ordered phase may have only one or several CDWs |l22| . 
For a rotational invariant system, it is often assumed that 
3 CDWs form a triangular lattice to minimize the break- 
ing of th e rot ational symmetry, as the 2D Wigner crys- 
tal state jl23l |. For systems with a strong lattice poten- 
tial, the system is often assumed to become an electron 
crystal state which preserves the point group rotational 
symmetry o f the background lattice, e.g. the rare-earth 
tritellurides[ll7|. 



For isotropic systems, outside the nematic phase, the 
p 3 term in Eq. (|4.9p f avor s that three CDWs form by 
a first-order transition [122]. However, inside a nematic 
phase, as we will show below, the nematic order parame- 
ter, which is coupled to p 2 , favors only one CDW and will 
compete with the p 3 term. For a continuous quantum 
phase transition, p 3 will be a subleading perturbation 
compared to p 2 , at least close enough to the transition. 
Hence, the smectic phase, a unidirectional CDW, will be 
energetically favorable. On the other hand, in the case of 
a first-order transition, depending on microscopic details 
either the smectic phase or the state with three CDWs 
would be preferred. We represent these different possi- 
bilities in the schematic phase diagram shown in Fig. [TJ 

On a square lattice, due to the point group symme- 
try of the lattice, the electron crystal phase usually con- 
sists of 2 CDWs perpendicular to each other. The phase 
transition between this phase and the FL may be second 
order due to the absence of the cubic term p 3 which, in 
contrast to isotropic systems, is prohibited by momen- 
tum conservation. We have confirmed this structure of 
the phase diagram in a microscopic mean-field calcula- 
tion. However, at the multi-critical point where both the 
CDW modes and the nematic mode are critical (the (0, 0) 
point in Fig. [T]), the coupling between CDWs and the 
nematic order parameter (Eq. (|4.12p ) is relevant. This 
suggests a fluctuation driven first-order transition near 
the multi-critical point. Hence, this multi-critical point 
is essentially unreachable. 

In this paper, we study the nematic- smectic phase 
transition and the smectic phase using a weak coupling 
approach by perturbing about a FL state. This ap- 
proach is consistent provided the nematic phase is narrow 
enough in coupling constant space so that the nematic- 
smectic transition is not too far from the FL phase. 

It is useful to compare to the classical version of 
this problem. The theory of classical (thermal) melt- 
ing in two dimension s, the Ko sterlitz-Thouless-Halperin- 
Nelson- Young theory[124j,|l25| (see Ref.[5[), is a theory of 
a phase transition driven by the proliferation of topologi- 
cal defects: a dislocation unbinding transition in the case 
of melting of a 2D Wigner crystal (a triangular lattice) 
into a hexatic phase, and disclination unbinding tran- 
sition in the hcxatic-isotropic phase transition. (The 
case of the square lattice was discussed only recently 
in Ref. [126| |). The reason for the success of the classi- 
cal theory of melting in two dimensions is th at, as in 
all Kosterlitz-Thouless phase transitions [H Il27j . at finite 
temperatures the classical ordered state with a sponta- 
neously broken continuous symmetry is not possible in 
two dimensions. Instead, there is a line (or region) of 
classical critical behavior with exactly marginal opera- 
tors. The defect-unbinding phase transition appears as 
an irrelevant operator becoming marginally relevant. 

In the case of the quantum phase transitions in two di- 
mensions that we are interested in, there are no such ex- 
act marginal operators available at zero temperature, and 
hence, no lines of fixed points available. Thus, the T = 
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FIG. 1: (online color) Schematic phase diagram at T — as 
a function of Ajv and Aodw defined in Eqs. (|4.6 [l and (|4,ll| l. 
The cross point of the two dash lines is the multi-critical point 
Ajv = Acdw = 0. The red thick lines stand for first order 
phase boundaries. Other phase boundaries may be first or 
second order. More complex electron crystal phases are pos- 
sible, for example, an anisotropic electron crystal phase where 
more than one CDWs and nematic coexist, but they are be- 
yond the discussion of this paper. 

quantum phase transition is not triggered by a defect- 
unbinding operator becoming marginal, but instead by 
making the coupling constant of an irrelevant operator 
large (as in standard continuous phase transitions, classi- 
cal or quantum). Instead, the quantum phase transition 
is closer to Landau-type (or, rather, Hertz-Millis like) 
description in that it is governed (as we will see) by a 
quantum-mechanical analog of the celebrated McMillan- 



deGennes theory for a nematic-smectic phase transition 
in classical liquid crystals in three dimensions [2, The 
approach that we will pursue here does not contain much 
of the physics of strong correlations as it begins with a 
state with well-defined fermionic quasiparticles. It also 
does not treat correctly the tendency of strongly corre- 
lated systems to exhibit inhomogeneous states and phase 
separation. The only way to account for this physics 
correctly is to use the opposite approach, a strong cou- 
pling theory of quantum melting of the crystal and stripe 
phases, as advocated in Ref.^j. So far, this theory only 
treats the physics deep inside a stripe phases, and the 
theory of their quantum melting to a nematic phase does 
not yet exist. Thus, although from a strong-coupling per- 
spective it would be highly desirable to have such a defect 
unbinding theory of this quantum phase transit ion ( such 
a description does exist for an insulating system [128| but 
its extension to a metallic state is not available and it is 
highly non-trivial), we will pursue instead a Hertz-Millis 
approach [47l . liH . [49[ to this quantum phase transition. 



D. The electronic nematic-smectic transition 

Nematic order will remove the degeneracy of CDW 
modes in different directions and select one CDW. As 
a result, the Brazovskii CDW multi-critical point be- 
comes just a critical point. For simplicity, we assume that 
the nematic order parameter is small enough so that a 
Landau- type expansion still makes sense, which is equiv- 
alent to assuming that the system is still "close enough" 
to the nematic-isotropic QCP. However, as we will show 
later, the critical theory we get using these assumptions 
has the only form allowed by symmetry, assuming ana- 
lyticity. 

By symmetry, the coupling between the CDW and the 
nematic field is 



^ int = ~- 9 JWfJ S iV ®^e- 2ie V(fc-?^-^M-fc,-0)+h.c. (4.12) 

I 



whose tensor form is shown in Appendix [X] Here, 9k 
is the polar angle of k. This term is irrelevant in the 
isotropic phase, but in the nematic phase, where N gets 
the expectation value N; this term will be of the same 
order as 52, which was defined in Eq. (|4.10[) , and hence 
it becomes important. 

Inside the nematic phase the amplitude fluctuations of 
the nematic order parameter are gapped while the orien- 
tational fluctuations, the nematic Goldstone modes, are 
gapless, at least strictly in the absence of a lattice and 



other orientational symmetry breaking couplings. Thus, 
deep enough in the nematic phase it is possible to inte- 
grate out the gapped nematic amplitude fluctuations and 
derive an effective theory involving the gapless nematic 
Goldstone mode. However, as the nematic-smectic phase 
transition is approached, the gap of the fluctuations of 
the smectic order parameter will get smaller and will ap- 
proach zero at the QCP. Thus, in this regime, the nematic 
phase has low-energy "fluctuating stripes" . This regime 
is the analog of that in conventional liquid crystals where 
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the McMillan-deGennes classical theory applies . We will 
now see how this theory arises in the quantum case. 
The leading term in £ int of Eq. (|4.12|) will be 



■111 I j^Ncos(20 k )p{k,uj)p(-k, 



(4.13) 



This term will stabilize the density wave in cither x or 
y direction and destabilize the other, depending on the 
sign of g. As a result, the nematic order will select a 
special direction along which only one CDW will form. 
Past this phase transition the system will be in a smectic 
state, a unidirectional CDW. For simplicity, we assume 
g > 0, which selects Qs in the y direction. 

Only the density fluctuations close to k = ±Qs matter 
for the low-energy theory. We define a complex field <J>, 
describing the density fluctuations around k = ±Qs as 



= p{q + Q s ,uj), 



(4.14) 



where q is small. The real part of $(r, t) measures the 
density fluctuations. 

Under a spatial inversion, <&(r, t) will become its con- 
jugate field r, t). Hence, the term — z<&^<9t<£> + h.c. is 
not allowed in the Lagrangian, and the dynamical term 
for $ is at least quadratic in time derivatives. 

The cubic term of the field in the isotropic-CDW tran- 
sition vanishes in the nematic-smectic transition, due 
to momentum conservation. By expanding Eq. (|4.12[) 
around q ~ 0, k ~ Qs and 4>n ~ 0, we obtain 



S =S$ N + J drdt(\d t <S>\ 2 - C v \dy<S>\ 



- C x \(d x - i^-M®] 2 - A s |$| 2 - usffl 

(4.15) 

Here S<p N is the action of the nematic Goldstone mode 

defined in Eq. gU). 

The action of Eq. P~T5)) is just a 2D version McMillan- 
de Gennes theory of the nematic-smectic transition in the 
classical liquid crystals but with z = 1 quantum dynam- 
ics. The constants in Eq. (|4.15p are 



C x 



AgN 



C y = C 



A s = A CDW - 2gN, 



4g 2 



us = «cdw - ~t • (4.16) 

An 



Here As is the energy gap of $ field, which mainly comes 
from the CDW gap defined in Eq. (|4.1ip . The correction 
term —2gN comes from the nematic ordering. The us 
term comes from the interactions between CDWs and it 
gets a correction from the amplitude fluctuations of the 
nematic order, which has been integrated out. The ne- 
matic Goldstone field 4>m couples to the CDW field $ as a 
gauge field with a "charge" Qs/2. Here the two in the de- 
nominator comes from the fact that the nematic order pa- 
rameter has an angular momentum 1 = 2. This gauge-like 



coupling is required by the rotational symmetry since, 
under spatial rotation by a small angle 8, the fields trans- 
form as $ — > exp(iQs x 9) <I> and <Pn — > 4>n + 10 (for the 
angular momentum channel I). In fact, with the sym- 
metry constrain and the assumption of analyticity, the 
action we show in Eq. (|4.15p is the only allowed form for 
the effective low-energy theor y, pro vided the topological 
excitations of 4>n are ignored [l29(. Therefore, although 
we only keep linear terms of N in our calculations above, 
which is valid close to the normal-nematic critical point, 
the action in Eq. (|4.15[) will have the same form even 
deep inside the nematic phase. 

The theory with the effective action given in Eq. (|4. 1 5[) 
has a critical field <& and gapless Goldstone boson </>at. A 
naive mean-field theory would suggest that this is a con- 
tinuous phase transition. In the case of the theory of clas- 
sical liquid crystals, where the same naive argument also 
holds, Halperin, Lubensky and Ma [fOJ used the 4 — e ex- 
pansion to show that there is a run-away behavior in the 
renormalization group flows, similar to that of supercon- 
ducting transition coupled to a fluctuating electromag- 
netic field. They concluded that in both cases the transi- 
tion is probably weakly first order, a fluctuation-induced 
first-order transition. In other terms, in the classical the- 
ory the coupling of the smectic to the nematic Goldstone 
mode (which has the same form as a coupling to a gauge 
field) is relevant. To ascertain what happens in the case 
of the metallic nematic-smectic QCP we will also need 
to take into account the effects of the fermionic degrees 
of freedom. We will see that the fermionic fluctuations 
change the critical behavior in an essential way. 



E. The electronic smectic phase: a unidirectional 
CDW 

In the smectic phase, the amplitude fluctuations of the 
order parameter, (5$, are gapped but the phase fluctua- 
tions, are gapless, as required by the Ward-identity. 
This happens in systems for which lattice effects can be 
neglected, and hence are described formally in a contin- 
uum, or if the smectic order is sufficiently incommensu- 
rate. Therefore, upon integrating out the gapped ampli- 
tude fluctuations 5$, the effective low-energy theory of 
the Goldstone mode becomes 

= Jdrdt [ KQ (dtM 2 - m{dl^) 2 - K 2 (d v ^) 2 ] , 

(4.17) 

When we are close to the nematic-smectic critical point, 
the coefficients of this effective action are 

Ko = |$| 2 , «i = ^S^-, K 2 = C y \®\ 2 , (4.18) 
Ws 

with i> being the expectation value of the CDW order 
paramete r. The vanishing of the stiffness {d x 4>$) 2 term 
[l3Ct Il3l| is required by the Ward identity of rotational 
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invariance. Thus, an underlying lattice, which will break 
the continuous rotational symmetry down to its discrete 
point group, will lead to a non- vanishing stiffness. Nev- 
ertheless, in many cases and particularly away from sit- 
uations in which the FS is strongly nested, the break- 
ing of rotational invariance can be parametrically small 
enough that at low temperatures its effects to a first ap- 
proximation can be neglected and treated perturbatively 
afterward. 

A simple scaling analysis of the effective action of Eq. 
(|4.17[) shows that, at the tree level, the scaling dimen- 
sions of space and time [x], [y] and [t], are —1,-2 and 
—2, respectively. Although the time direction and the y 
direction scales in the same way, the x and y directions 
now scale differently. This is a typical phenomenon for 
anisotropic states. Now the effective dimensions of this 
quantum theory is 5. Hence, our theory is above its (up- 
per) critical dimension. So the higher-order interactions 
of the smectic Goldstone mode will be irrelevant, if we 
don't consider topological defects. The fact that we are 
above the critical dimension also tells us that the quan- 
tum fluctuations in the quantum smectic phase will not 
destroy the long-range order. 

This scaling is very different from the classical smec- 
tic phase of 3D, where [x] = [y] = — 1 and [z] = —2, 
if the modulation is on the z direction. This classical 
theory is at its lower critical dim ension, an d long-range 
order is destroyed by fluctuations [l30l Il3l| , resulting in 
a power-law quasi-long-range order. This system has a 
line of critical points, so the higher order terms of the 
action that need to be considered were found to l ead t o 
logarithmic corrections to the power-law behavior jl32| . 

The above analysis implies that our quantum prob- 
lem is above the lower critical dimension. Therefore all 
these effects of the 3D classical smectic phase will not 
be present in the 2D quantum case. The scaling behav- 
ior of a 2D quantum system is similar to the columnar 
state of the classical liquid crystals, instead of that of 
classical smectics. The classical columnar state has two 
density waves so that it is a solid in two directions but a 
liquid in the third direction. The Goldstone fluctuations 
of this state scale as [x] = —1 and [y] = [2] = -2 |, 
which is the same as in the present case, if we consider 
the time direction in our problem as the z direction. The 
difference between the classical columnar state and the 
2D quantum smectic state is that in the 3D columnar 
state, the Goldstone mode is a planar vector but in the 
present case it is a scalar. 



V. COUPLING THE ORDER PARAMETER 
THEORY TO FERMIONS 

We will now proceed to couple the phenomenological 
theory of the nematic and smectic phases to a system of a 
priori well-defined fermionic quasiparticles described by 
the Landau theory of the FL. In a fermionic liquid crystal 
state, the bosonic order-parameter fields, defined above, 




(c) Q S < 2k F (d) Q S = 2k F 



FIG. 2: (online color) The FS of the nematic phase (a and 
b) and the reconstructed FS in the smectic phase (c and d). 
(a) and (c) are for Qs < 2k F at the QCP and in the smectic 
phase respectively, while (b) and (d) are Qs = 2k F , also at 
the QCP and in the smectic phase respectively. In (a) and (b) , 
the black dots marked the non-FL points on the FS caused by 
the smectic mode fluctuations at the nematic-smectic QCP. 
The relevance of the points in (c) is explained in Sec. IVIII 
In (c) we have show the case of Qs to be comparable to 2kF 
so as to keep the FS reconstruction simple. Here we show 
the effective Brillouin Zone with an open orbit and a closed 
pocket. The reconstructed FS of case (d) is partially gapped 
and the FS has an open orbit. 



will couple to the fermions. 

Let us define i/j' (x, t) and ip(x, t) to be the fermion cre- 
ation and annihilation operators of a FL. We will assume 
that the FL has a well-defined FS, which for simplicity 
we will assume is circular. (For lattice systems the FS 
will have the symmetry of the point group of the lat- 
tice.) The Fermi wave vector is kp. The Fermi velocity 
is set to 1 so that the energy and momentum have the 
same units. Consisten t wit h the assumptions of the Lan- 
dau theory of the FL |l33l ] the effective Hamiltonian of 
the fermionic quasiparticles will be taken to be that of 
a free Fermi system, with a well-defined FS, and a set 
of quasiparticle interactions parametrized by the Landau 
parameters. These interactions are irrelevant in the low- 
energy limit of the FL but play an important role in the 
physics of electronic liquid crystal phases (33[. In any 
case in our discussion it will be unnecessary to include 
the Landau parameters explicitly since their effects will 
already be taken into account through the coupling to 
the liquid crystal order parameters. 

By symmetry, the nematic order-parameter field, N, 
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couples to the fermion density quadrupole [3c 

Qm = ^(r,t) ( %Jy _f% 2 ) mt). 



(5.1) 



In 2D, since the rotational group is 5*0(2), the density 
quadrupole can be defined in terms of a two-component 
real director field (i.e. a headless vector) or, in terms of 
complex field 



Q(r,t) = Qn(r,t)+iQ 12 (r,t) 

= ^m {dx+ J dy)2 mt). (5.2) 

Same as the nematic order parameter, Q is also invariant 
under rotations by 7r. 

The coupling between Q and N is 



g N / drdt(Q*N + h.c). 



(5.3) 



Here g^ is a coupling constant. Again, the chiral sym- 
metry of the system requires that the effective action 
depends only on the real part of Q^N, and that there is 
no dependence on the imaginary part, since it is a pseudo 
scalar. The tensor form of this coupling is shown in Ap- 
pendix In what follows we choose the sign of gjy to be 
negative, so that a positive expectation value of the ne- 
matic order parameter TV means a FS stretched along the 
x direction and compressed in the y direction, as shown 
in Figs. E^a) and (b). 

The sign of gjv alone is not important. What matters 
is the relative sign between gjy and the coupling constant 
g defined in Eq. (|4.12p . Under a redefinition of N be- 
coming — N, both gpf and g change sign. If g x g^ > 0, 
Qs prefers the direction in which the FS is stretched, but 
when g x g^ < 0, it prefers the direction where the FS is 
compressed. In general, the sign of g x g^ is determined 
by microscopic details of the system to which this model 
may apply. 

If Qs = 2fcp, very close to a nesting condition the 
curvature of the FS controls the CDW instability as it 
controls how singular the charge susceptibility is near 
the nesting wave vector. In this case one finds that it 
leads to the condition g x g^ < 0, when Qs connects 
two points on the Fermi surface where the curvature is 
smallest, as shown in Figf2Jb). In general, far from a 
nesting condition, the curvature of the FS alone is not the 
dominant factor, and the sign of g x g^ may be positive 
or negative, depending on the microscopic details. 

The smectic order-parameter field should be coupled 
to the CDW of the fermions. The CDW operator of the 
fermions, close to the ordering wave vector Qs, is 



n(q, uj) 



AhAO 

3 1p (k + Q S + q,n + to) $(k, fi), 



(2tt) 2 



where q ~ 0. The smectic order-parameter field $ cou- 
ples to this fermion density wave n as 



gs / drdt(n>$ + h.c.) 



(5.5) 



Integrating out the bosons, attractive four-fermion inter- 
actions are generated of the form 



#M 2 + #H 2 . 

Aw As 



(5.6) 



Hence, the order-parameter fields can be regarded 
as Hubbard- Stratonovich fields used to decouple four- 
fermion interactions. In this picture, the couplings be- 
tween the order-parameter fields and fermions are mea- 
suring the strength of the attractive four-fermion term. 

Gapless fermions will introduce nonanalytic terms to 
the low-energy effective theory of the nematics and 
smectics. For the case of the nematic order param- 
eter, it was shown by Oganesyan and co-workers [33[ 
that the fermions generate nonanalytic Landau damp- 
ing terms [13, |48j], so the theory of the isotropic-nematic 
metallic QCP becomes 



S N - g 2 N N(0) 



dqduj 



-K N q 2 )N^(q,uj)N(q,Uj), 
(5.7) 



The nematic susceptibility at this FL-nematic QCP is [331] 
X N (q,u)=-i(N^q,u)N(q,u))r e t 

(5.8) 



.'/v -Villi \ — - /• \ <f 



The phase mode of the nematic order-parameter field 
in the nematic phase, the nematic Goldstone mode, has 
an effective action of the form 



5V N = g 2 N N 2 N(0) 



x / ^ (^sin 2 2 lPq -K 1 q 2 x -K 2 q 2 y ) |<M<Z»| 2 

(5.9) 



(2tt) 3 



where N is the expectation value of the nematic order 
parameter and (p q is the angle between q and the main 
axis direction of the nematic ordering. The stiffnesses 
K\ and K 2 (the Frank constants) are given in Ref.[33(. 
With this action, it follows that the transverse nematic 
susceptibility in the electron nematic phase is[33| 

XxOf> w ) = - iN 2 {(j> N (q,u))(j) N (-q,-uj))^ t 



g 2 N N(0) ( j sin 2 2 Vq - K lQ 2 x - K 2 q 2 y 



(5.4) 



(5.10) 
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For the case of a nematic order parameter aligned 
along the x-axis, the angular factor becomes sin 2 2ip q = 

For the case of a charged smectic, a unidirectional 
CDW, a similar effect will be observed. Besides, if Qs 
connects to points on the FS which have just the opposite 
Fermi velocity as shown in Fig. (Hb), the discontinuity 
leads to another type of nonanalytic terms as will be 
shown in Sec. IVI Bl 



VI. THE NEMATIC-SMECTIC METALLIC 
QUANTUM CRITICAL POINT 

In this section, we study the metallic nematic-smectic 
QCP. The two cases shown in Figs. HJa) Qs < 2kp, and 
(b) Qs = 2fci?, are studied separately. 

Deep in the nematic phase, the amplitude fluctuations 
of the nematic order parameter are gapped, and the low- 
energy fluctuations are due to the nematic Goldstone 
mode, (f>N, whose action is given in Eq. (|5.9p . How- 
ever, as the nematic-smectic QCP is approached (from 
the nematic side) the fluctuations of the smectic order 
parameter become progressively softer and, provided the 
quantum phase transition is continuous, become gapless 
at this QCP. In this scenario, the nematic phase looks 
like a "fluctuating stripe" phase qualitatively similar to 
the phenomenology of the cuprate superconductors, as 
discussed in Sec. IIII Al 

The case Qs > 2kp will not be discussed here. The 
reason is that since now the CDW fluctuations with 
Qs > 2kp cannot decay into particle- hole pairs, in 
this case fermions only renormalize the coefficients of 
the smectic effective action, while the nematic fluctua- 
tions will still be Landau damped. For isotropic sys- 
tems and for Qs > 2kp, the CDW (Lindhard) suscep- 
tibility x(Qs) m general decreases faster than linear as 
x(2kp) — A \/Qs — 2kp, where A is a constant. This im- 
plies that a CDW with Qs > 2k f is unlikely to be realized 
as it would require an anomalously attractive interaction 
at a large Qs- However, for a lattice system the phase 
fluctuations of the nematic mode 4>n get gapped by lat- 
tice anisotropics and in this case the fermions only yield 
the trivial effect of renormalizing the coefficients of the 
effective action at the CDW transition. 



A. Q s < 2k F 

For Qs < 2kp, the leading contribution to the effec- 
tive action of the order parameter field, resulting from 



integrating out the fermions, has the form 
dqduj 



9% 



(2n)' 



U(Q s + q,^) |<%»|' 



(6.1) 



Here H(Qs + q, ui) is the CDW susceptibility of the 
fermions, given by the fermion loop integral (bubble) 



IL(k,oj) 



dp 



np[e(p + k)] - n F [e{p)} 



lu — e [p+ k) + e {p) + i0 + sign(w) 

(6.2) 



where np(k) is the Fermi-Dirac distribution function. 

The static part of the fermion CDW susceptibility de- 
pends on the details of the dispersion relation from way 
above the FS to the bottom of the band. However, since 
n(fc, oj = 0) is analytic for k < 2kp, the static part, 
H(Qs + q,uj = 0), will not change the analytic struc- 
ture of Eq. (|4.15p . but just renormalize the coefficients, 
in particular the critical value of the coupling constant. 
The important contribution comes from the dynamical 
part, n(<2s + q, uj) — H(Qs + q, 0). The singular contri- 
butions to this integral are dominated by the behavior of 
the integrand around the four points on the FS, which 
are connected by the ordering wave vector Qs, as marked 
with black dots on Fig. [2ja). If we expand the disper- 
sion relation of the fermions around these four points, 
e(q) = ±v x q x ± v y q y , to leading order we get a Landau 
damping contribution 



U(q + Q s ,uj)-Il^+Qs,0) 



2-!TV x Vy ' 



(6.3) 



which is linear in |cj|. The formula above can be checked 
by taking the limit of Qs <C 2kp or Qs < 2kp. In 
these two regimes, the fermion loop integral can be com- 
puted by RPA without expanding the dispersion rela- 
tions around the four points. After setting vp = 1, 
for Q s < 2k F , one finds U(q + Q s ,u) = iN(0)\uj\/Q s 
with N(0) being the density of states and for Qs < 2kp, 
n(<f + Qs,u) — i\u}\y/kp / (2ny / 2kp — Qs), which can be 
reached by expanding Eq. (|6.9p . Both of them agree 
with the general formula given above. 

The term linear in u in the effective action for the 
smectic field $ of Eq. (|6.3p . which is due to the con- 
tributions of the fermions, dominates over the "naive" 
dynamical term proportional to u> 2 of the phenomeno- 
logical theory. We can thus write an effective action for 
the electron nematic-smectic quantum phase transition 
of the form 
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<iqdUJ C i\co\ \<S>{q,u)\ 2 - J drdt (^Cyldy^ + Cxlidx-i^^ + Asl^+usm* 



(2tt) 3 
dqduj 
(2^)3 



i^o— sin 2 2<p q - K x q 2 x - K 2 q 2 y \<j) N (q,ij)\ 



(6.4) 



where C = g 2 /{2nv x v y ), K = g 2 N N 2 N(0) and K lf2 = 
g 2 N N 2 N(Q)Ki a - The point A s = and u s > is the 
nematic-smectic critical point. With the nonanalytic dy- 
namical term, the dynamic critical exponent of the field 
4> becomes z = 2, instead of z = 1 as it would generally 
be in the absence of fermions (or, if the fermions were 
gapped as in the case of an insulator). 

The nematic Goldstone mode 4>n has a dynamic criti- 
cal exponent z = 3 [Hj], larger than the z = 2 exponent 
for the smectic fluctuations. Thus, the Goldstone mode 
of the nematic order parameter cf>N and the smectic $ 
fluctuate on very different energy scales, with 4>n being 
the low-energy mode. If we only focus on the asymp- 
totic low-energy theory, we should integrate out the high- 
energy mode $. This process will lead to an effective 
theory of <pN ■ In turn, the low-energy mode 4>n will me- 
diate interactions of the field <!>. However, we will show 
by a scaling argument that in the case of the quantum 
metallic system the coupling between the smectic field 
and the nematic Goldstone mode is irrelevant. 

The action of Eq. (|6.4p is invariant under a rescaling 
parametrized by a factor b 



t -> b~ 3 t, 

$(r,*) -► & 3/2 $(&-y,k- 3 f), 



c 



x,y, 



c 

A s 6 2 As, 

#0,1,2 - b 3 Kn 



r — > b f, 
Co — > b 1 Co, 
Qs -> bQ s , 
us -> b~ 1 u s , 



(6.5) 



where C x , C y , Kq, K±, and K% are the stiffness in Eq. 
(EH). 

When As = 0, at the tree level and in the long- 
wavelength regime, both the gauge-like "coupling con- 
stant" Qs and #0,1,2 scale to infinite, but the ratio 
Ql/ifo.1,2 scales to as a function of i> -1 . This im- 
plies that the gauge- like coupling is irrelevant. Quantum 
fluctuations may change the tree-level scaling behavior as 
we include loop corrections. However, for large enough 
#0.1,2 or small enough Qs-, the irrelevancy of the gauge- 
like coupling will not be changed. As a byproduct, we no- 
tice that Co and us scale to zero in the long-wavelength 
regime, which means that these two terms are irrelevant 
also. However, we should keep in mind that these oper- 
ators are actually dangerous irrelevant, in the sense that 
Co is necessary to find the proper equal-time correlation 
function for $ and us is necessary for stability in the 
ordered phase, and they are only irrelevant at this QCP. 



Notice that at the QCP there are two critical modes: 
the amplitude of the CDW order parameter, which has 
z = 2, and the transverse (Goldstone) mode of the ne- 
matic phase, which has z — 3 (and it is clearly dominant 
at low enough energies). Thus, we also need to check 
the scaling behavior of t — » b~ 2 t and r — > b^t for the 
high-energy mode. Under this rescaling, 



$(r,i) -> 6$(6"V,6~ 2 t), 
Qs -> bQ s , 
us -> u s , 

#1.2 -» & 2 #1.2. 



As - 6 2 As, 
#0 -» b 3 K , 

(6.6) 



At the critical point where As = 0, it can be seen that 
Q 2 s /Kq also scales to as in the long-wavelength 
limit, which means, for the z = 2 mode, the gauge-like 
coupling is still irrelevant. 

These conclusions are confirmed by one-loop perturba- 
tion theory calculations, presented in Appendix [Bl where 
we show that integrating out $ (or does not change 
the action of 4>n (or <&). This is one of our main results. 

In conclusion, there are two essentially decoupled soft 
modes at the nematic-smectic QCP. The nematic Gold- 
stone mode is governed by the same action as in the ne- 
matic phase, Eq. (|5.9p . Since as the nematic-smectic 
QCP is approached from the nematic side, the nematic 
Goldstone mode and the smectic order parameters effec- 
tively decouple; the effective action for the smectic field 
$ in this limit reduces to 



Ss 



dqduj 



(iCoM - C x q 2 x - C y ql - As) |4>(g»| 2 , 

(6.7) 



which implies that the dynamic smectic susceptibility is 

(6.8) 



X 



1 



iC U - C x q 2 - C y q 2 - As 



where q is the momentum measured from the ordering 
wave vector Qs- Here As > on the nematic side of 
this QCP. On the smectic (stripe) side of the quantum 
phase transition the $ 4 coupling, which is (dangerous) 
irrelevant at this QCP, cannot be ignored as it stabilizes 
the smectic ground state. The smectic susceptibility in 
the ordered smectic phase differs from that of Eq. (|6.8p in 
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FIG. 3: (color online) The spectral density of the smectic sus- 
ceptibility at the nematic-smectic QCP, Imxs, as a function 
of q and to for Qs < 2fci?. The spectral density is singular 
near the origin (lower right corner) and decays monotonically 
away from there. Here we show contour plots at constant 
spectral density with values from up to 2000. The red line, 
uj — Cxq't/Co + C y qy/Co, marks the peak of the spectral 
density as a function of momentum q parallel to the nematic 
orientation. The inset is the energy dependence of Imxs at a 
fixed small momentum (along the dashed vertical line). 



two standard ways: a) it acquires the usual delta function 
term peaked at the ordering wave vector, & 2 5(q), where 
<f> is the expectation value of the smectic order parameter 
and q is measured from the ordering wave vector Qs, and 
b) the "mass term" (As) in the denominator of the sus- 
ceptibility becomes 2|As|. At the nematic-smectic QCP, 
As = 0, the smectic fluctuations are described by an 
overdamped critical mode with z = 2. The spectral den- 
sity of xs at the QCP is shown in Fig. [3] It shows that 
most of the spectral weight is at small lo and q, which is 
typical for a critical mode, and the energy distribution 
curve at fixed momentum has a broad peak marked by 
the red line (lo = C x q 2 /Co + C y qy/Co), indicating an 
overdamped critical mode with z = 2. 

The low-energy physics of the system will be domi- 
nated by cj>Mi the Goldstone mode of the nematic field, 
whose behavior has been studied extensively in Ref. 
[HI, On the other hand, at higher energy scales 

(u; ~ q 2 ), the effects of field the amplitude mode of the 
smectic (or stripe) fluctuations, will become observable 
and, in this range, the system effectively has "fluctuating 
stripes" . Notice that if the nematic Goldstone mode be- 
comes gapped, say by the effects of the lattice, the smec- 
tic amplitude fluctuations become the only low-energy 
modes left. Finally, since the field <J> at the critical point 



has 2 = 2, the effective total dimension of this theory is 4. 
A standard Hertz-Millis type argument [4?l \I§L |49| implies 
that in this case the |$| 4 term is marginally irrelevant. 
Therefore, the Gaussian fixed point will have the cor- 
rect scalings, up to logarithmic corrections. In particular, 
contrary to what happens in the classical case where this 
transition becomes weakly first order, at the nematic- 
smectic QCP the electron smectic susceptibility will ac- 
quire logarithmic corrections to scaling in the ui — > and 
q — > limit, as can be deduced from standard ar gume nts 
in classical and quantum critical phenomena [49lTl34l | . 

The gapless smectic amplitude mode $, the fluctuating 
stripe, can be detected in inelastic light scattering exper- 
iments, much in the same way as in the case of conven- 
tional CDW materials. The existence of this mode also 
has observable effects on thermodynamic properties such 
as the low-temperature heat capacity. Since the dynamic 
exponent now is z — 2, the fluctuation of the amplitude 
of the smectic mode leads to a C ~ T contribution, which 
is subleading compared to the C ~ T 2 / 3 contribution of 
the nematic Goldstone mode 4>n- 

Besides the subleading contribution to the heat ca- 
pacity, the critical fluctuations of the smectic order pa- 
rameter also have an observable effect on the fermions 
as they profoundly change the character of these exci- 
tations. As pointed out by Oganesyan and coworkers 
[HI, the fermionic states in the nematic phase will be- 
come a non-FL due to the effects of the fluctuations of 
the nematic Goldstone mode, an overdamped mode with 
dynamic critical exponent z = 3. To leading order in per- 
turbation theory, these authors found that for the most 
of the FS, the imaginary part of the fermion self-energy 
(the quasiparticle scattering rate) T,"(kp,Lo) ~ |w| 2 / 3 , 
which, as lo — ► (i.e. as the FS is approached) vanishes 
slower than lo. Thus, in almost all of the FS the fermionic 
quasiparticle is no longer a well-defined state as the quasi- 
particle pole in the fermion Green function is lost. This 
is the signature of a non-FL. However, for fermionic ex- 
citations propagating along the four main-axis directions 
of the nematic FS, the quasiparticle scattering rate now 
scales as T/'(kp,uS) ~ |w| 3 / 2 . Although this is not the 
conventional lo 2 behavior expected in an ordinary FL, 
nevertheless it is still consistent with the existence of a 
pole on the fermion spectral function, and a well-defined 
quasiparticle. Hence, in the nematic phase, except along 
these four special directions, the quasiparticles are not 
well defined. 

The |w| 2 / 3 behavior of the quasiparticle rate is in clear 
conflict with conventional FL behavior. It also implies 
that perturbation theory is breaking down in this sys- 
tem. Two approaches have been proposed to assess the 
non-perturbative behavior of the system. Using the non- 
pert urba t ive approach of h igher dimensional bosoniza- 
tion [HE EM UM EH, EH, Lawler and Fradkin showed 
that at the non-perturbative level the \lo\ 2 ^ 3 perturbative 
correction to the quasiparticle rate leads to a dramatic 
change in the behavior of the fermion propagator, which 
they found to have a vanishing quasiparticle residue and 
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to exhibit a form of "local quantum criticality" as it 
scales in frequency but not in momentum (45l. [4(1 . On 
the other hand, Chubukov and Khveshchenko jl40l | used 
a resummed perturbation theory approach (on a sim- 
ilar problem) and argue that the |w| 2 / 3 behavior per- 
sists to all orders in perturbation theory. Although it 
is presently an open problem how to reconcile these two 
results, both analysis lead to the conclusion that the 
fermionic quasiparticles do not exist as well-defined ex- 
citations at the nematic-FL QCP and throughout the 
nematic phase (provided the nematic Goldstone modes 
remain gapless and overdamped). 

At the nematic-smectic QCP, the fluctuations of the 
smectic field $ (i.e. the "fluctuating stripe" mode) will 
also contribute to the fcrmion self-energy corrections (Sec 
Appendix IE II for details). For most points on the FS, the 
contribution to the quasiparticle rate H"(kF,oj) of the 
fluctuations of the smectic field $ is proportional to oj 2 , 
which is consistent with a conventional Landau behav- 
ior and a well-defined quasiparticle. However, at special 
points on the FS satisfying e(fc) = e(k + Qs) (shown 
in Fig. (Ha) as the four black dots), the contributions 
of the fluctuations of the smectic order-parameter field 
$ to the quasiparticle rate scale as Y^"(kp,uj) ~ M 1 / 2 . 
This M 1 / 2 behavior dominates even over the |w| 2 ^ 3 con- 
tribution of the nematic Goldstone mode. Hence, at the 
nematic-smectic QCP the quasiparticle residue will van- 
ish at these special points of the FS. Note that these spe- 
cial points are precisely the positions on the FS where FS 
reconstruction will take place due to the development of 
CDW order, as shown in Figs. [U (a) and (c). Hence, it 
is not surprising to see that strong derivations from the 
Landau FL picture appear at these points at the QCP. 
This non-FL behavior is just the prelude of the FS re- 
construction in the ordered phase. 

Finally, if the continuous rotational symmetry is bro- 
ken explicitly by the anisotropic effects of the underlying 
lattice (say through an anisotropic band structure) or by 
external fields, at very low energies the nematic Gold- 
stone mode of the nematic phase will have a finite gap 



(generally z = 1). If this gap is small enough, sufficiently 
close to the FS the fermion quasiparticle rate will show 
a crossover from the above mentioned |w| 2 / 3 above this 
gap to a conventional uj 2 characteristic of a Landau FL. 
(Note that at the nematic-FL QCP the M 2 / 3 behavior 
is still obtained even for a lattice system [36|.) Neverthe- 
less, at the nematic-smectic QCP, the quantum critical 
fluctuations of the smectic mode still generate M 1 / 2 cor- 
rections to the quasiparticle rate of the fermions at the 
special points where the FS is going to be reconstructed. 



B. Q s = 2k F 

We now consider the special case of Qs — 2kp. For 
the same reason as mentioned above, the coupling be- 
tween the nematic Goldstone mode and the smectic field 
is irrelevant. Hence, at sufficiently low energies and close 
enough to the QCP, we can ignore their coupling and con- 
sider the effective theory of the smectic order-parameter 
field alone. The effects of the irrelevant coupling to the 
nematic Goldstone mode can be put back in perturba- 
tively a posteriori. 

When Qs = 2kp, Eq. (|6.1|) is still valid, although the 
structure of the fermion loop integral is very different. 
We can compute the effective theory of the bosonic modes 
by evaluating the fermion loop integral in the same way 
as we did in the previous section for the Qs < 2kF case. 
Again, we find that Tl(Qs, 0) depends on the microscopic 
details of the fermion dispersion relation, but Tl(Qs + 
q, ui) — H(Qs,w = 0) are dominated by the behavior of 
the integrand around the two points on the FS connected 
by Qs as marked in Fig. [2Kb). We expand the dispersion 
relation around these two points as e = fi±5k y + Kdk 2 /2, 
where fj, is the chemical potential, k is the local curvature 
of the FS and 8k is the momentum measured from these 
two points. By evaluating the fermion loop integral using 
this dispersion relation, we can determine the low-energy 
Lagrangian density of the field $ to quadratic order 



U(q,w) = - [CUjQy + ^+U +)Jqy + ^-uj+a 1 q y + a 2 q 2 x j \<P(q,uj)\ 2 , (6.9) 

which implies that the smectic susceptibility at the QCP of this case is 

X s (g» = -i(&(q,cj)$(q,u)) Iet = . ■ ■ ~* =^ . (6.10) 



C, q y + ^ + u> + i0+ + \J 'q y + ^ - U! - i0+ j + a x q y + a 2 q 2 



r 



Here £ = g s /(2it^/K) and are determined by mi- 
croscopic details. Notice that ( diverges as k vanishes, 
which is related the fact that a flat FS has a logarith- 



mically divergent CDW susceptibility. In Sec. IVIDI we 
will show that at an inflection point of the FS, where 
K = 0, there is a stronger non-analytic behavior of the 
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form ^/q y + lo. 

For this Lagrangian density, a naive dimension count- 
ing suggests the scaling dimensions [q x ] = 1, [q y ] — [lo] = 
2. Under this scaling, a\ and cti are irrelevant. However, 
if we take a small lo expansion for q y + Kq y /4 < 0, the 
Lagrangian density becomes 




(6.11) 



since y q y + Kq%/4 + lo and y q y + nq x /4 — lo cancel each 

other in the u> — > limit when q y + Kq x /4 < 0. This 
Lagrangian density has scaling law [q x ] = 1, [q y ] = 2, 
and [lo] = 3. This scaling is only valid inside the 
particle- hole continuum, while the naive scaling, [q x ] = 1, 
[q y ] = [lo] — 2, is valid outside. This different behavior 
arises because Qs = 2fcp is located at the edge of the 
particle-hole continuum. Hence, the long wave length 
fluctuations may be inside or outside the particle-hole 
continuum, which leads to two possible different scaling 
behaviors. Among these two different scaling behaviors, 
the [q x ] = 1, [q y ] = 2, and [lo] = 3 is the low energy mode 
in the long-wavelength limit. Hence, this mode domi- 
nates the low-energy physics and the scaling behavior. 

The fact that q x , q y , and lo have different scaling di- 
mensions is typical of anisotropic systems. For instance, 
in the smectic (stripe) phase of the quantum Hall state, 
a similar scaling was found in Ref. [97|. ( This scaling 
behavior was later on proved to be unstable [l4l| due to 
the existence of an infinite set of marginal operators in 
that system.) Although the problem we are discussing 
here and the smectic quantum Hall state share the same 
scaling dimensions, they are actually quite different. In 
the case of the theory we are discussing in this section, it 
is the theory of the QCP of the nematic-smectic transi- 
tion, while in the quantum Hall case, the same scaling is 
found in the smectic phase. Second, the different scaling 
dimensions in the x and y directions in our case are due 
to the existence of the FS, which results in two different 
scaling dimensions depending on whether the momentum 
is perpendicular or parallel to the FS. In contrast, in the 
quantum Hall smectic phase, it is due to the residue sym- 
metry of the broken rotational symmetry, which is the 
same as in the smectic phase of the classical liquid crys- 
tals. Third, the scaling dimension 3 in the time direction 
is due to the nonanalytic dynamical term in the present 
problem we are discussing, but in the quantum Hall smec- 
tic, it is due to the explicitly broken time-reversal sym- 
metry. Most importantly, the quantum Hall smectic is 
essentially an insulator in the direction perpendicular to 
the stripes. In our case, however, the stripe has not yet 
formed at the QCP, and the system is an anisotropic con- 
ductor in all directions. For our system, even inside the 
smectic phase, the FS is just partially gapped, so the 
conductivity in the direction perpendicular to the stripes 



is still non-zero. This difference is very important, since 
it means that in the quantum Hall smectic, the system 
is actually a II? sliding Luttinger liquid, but in our case, 
the system is a full 2D structure. It is precisely the ex- 
istence of the sliding symmetry that makes the [q x ] = 1, 
[q y ] = 2, and [lo] = 3 scalings unstable for the quantum 
Hall smectic phase. For our system instead, due to the 
absence of the sliding symmetry, our [q x ] = 1, [q y ] = 2, 
and [lo] = 3 scalings will not experience the same insta- 
bility as the quantum Hall smectic phase. 



On the other had, non-Gaussian terms play an impor- 
tant role in our problem. The leading non-quadratic term 
is not the naive |$| 4 term as in an ordinary <i> 4 bosonic 
theory. Instead, the discontinuity at the FS introduces 
a non- analytic term |<I>| 5 / 2 , as shown in Appendix [Ul A 
similar non-analytical term has been found in the FL- 
ferromagnets transition by Maslov, Chubukov, and Saha 
[l42| . where a non-analytic term |m| 3 in the thermody- 
namic potential, with m being the ferromagnetic order 
parameter, results in non-analytic behaviors for the spe- 
cific heat coefficient and the spin susceptibility. It is also 
reported in the same reference that this \m\ 3 term will be 
replaced by a term cx |m| 7 / 2 at the QCP. Our theory will 
have effective dimension 6 at the Gaussian fixed point 
with scaling [q x ] — 1, [q y ] — 2, and [lo] = 3, where |$| 4 is 
irrelevant ([$] = 2), but |$| 5 / 2 is relevant. This relevant 
|$| 5 / 2 term changes the scaling behavior of the critical 
theory away from the Gaussian theory. 



Hence, in contrast with the case Qs < 2kp, the quan- 
tum critical theory at Qs = 2/ip is not controlled by 
the Gaussian fixed point. This conclusion agrees with 
the results of the "small momenta" regimes (Qs = 2/c_f 
incommensurate CDW critical point) of the theory of a 
CDW-FL QCP of Altshuler et al %2%, who also found 
a relevant perturbation at the Gaussian fixed point and 
hence a runaway RG flow. Therefore, the Gaussian the- 
ory fails. As a result, the irrelevant terms, which we 
ignored in this study, will determine the fate of this tran- 
sition. In general, there are two possible situations. If 
the runaway RG flow has a (non-perturbative) disconti- 
nuity fixed point, the phase transition will become first 
order. This scenario is known as the "fluctuation-driven 
first order transition". However, there is also the possi- 
bility that the runaway RG flow has a non-trivial (and 
also non-perturbative) quantum critical fixed point, in 
which case the transition is still second order but has a 
different scaling. It may even be possible to change from 
a first-order transition to a second one by tuning some 
control parameters and going through a quantum tricrit- 
ical point. In any case, although it is generally assumed 
that a fluctuation-induced first-order transition normally 
results, actually it is not possible to determine which one 
of these two scenarios actually holds by theperturbative 
arguments we are using here (and in Ref.[28|). 
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C. Commensurate CDW on a lattice of Ref. [28| , In this case, the Gaussian part of the effec- 

tive Lagrangian density will become 

Our analysis can be used for the case of a commensu- 
rate CDW QCP as well, i.e. the "large momenta" regime 



U = -C\^q y + ^f+uJ + ] l qv + ^-uJ+ \/-q y + + sj-q y + ?f - wj |<%»| 2 , (6.12) 

which leads to the smectic susceptibility 

-vc 



q y + ^ + w + i0+ + ^q y + ^ - u - i0+ + ^-q y + ^ + lu + i0+ + -q y + ^-^-10+ 

(6.13) 



The last two terms in the Lagrangian density appear, 
due to the mirror effect between momentum k and G — k 
introduced by the band structure, where G = 2Q$ is 
a reciprocal lattice vector. The static part of the La- 
grangian density docs not have a full cancellation as the 
incommensurate case in Eq. (|6.11[) due to the fact that 
±q y + Kq x /4 can not be negative at the same time. Hence, 
the naive scaling law, [q x ] = 1, [q y ] = 2, and [lu] = 2, 
is always valid, which makes the effective dimension of 
this theory 5. At the tree level-since $ has dimension 
2, |$| 5 / 2 is marginal. In Ref. [28(], it is shown that in 
this regime the coupling between the CDW order param- 
eter and the fermions is marginal at the tree level and, 
at one- loop level it becomes marginal irrelevant. There- 
fore, the |$| 5 / 2 term should be marginal irrelevant here 
too. We conclude that the quantum critical theory of 
the commensurate CDW critical point is Gaussian, with 
logarithmic corrections to scaling. 

In this case, the contribution of the quantum criti- 
cal smectic field to the low-temperature specific heat is 
proportional to T 3 / 2 . The quasiparticle scattering rate 
Y,"(kp,Lu) acquires a contribution (due to the quantum 
smectic fluctuations) throughout most of the FS propor- 
tional to lu 2 , which can be ignored in low-energies and 
it is consistent with a finite quasiparticle pole. However, 
for the fermionic excitations close to ±kpe y , '(hp, lu) ~ 
\lo\, which suggests a marginal FL behavior. Details of 
this analysis can be found in Appendix IE II The result 
Ti" (kp,u>) y \lu\ for k = ±kpe y agrees with the findings 
of Ref. [2|. ' 



D. The special case of a FS with inflection points 

In the case of a two-dimensional system the FS is a 
curve. The point on the FS where the curvature van- 




(a) (b) 



FIG. 4: (online color) The FS near the inflection points, (a) a 
FS with 4 inflection points and only the two connected by Qs 
will be considered; (b) the FS near the two inflection points 
by shifting them together and rotated. 

ishes is geometrically an inflection point (Fig. HJ. In 
the Qs = 2kp case, if the two points connected by Qs 
happen to be two inflection points, the scaling behavior 
will become even more non-analytic than the cases we 
discussed above. Very close to the FS, we can assume 
that near the inflection points the dispersion relation of 
the fermions has an expansion of the form 

e{q + Qs/2) = fi + q y + aq 3 x + bq 4 x + . . . (6.14) 

Here /i is the chemical potential, the wave vector q has 
components q x and q y , which are the components per- 
pendicular and tangent to the FS, all measured from the 
inflection point, and a and b are two constants. Assum- 
ing that the band structure has inversion symmetry, to 
obtain the dispersion relation around the other inflection 
point, we just need to reverse the vector q. The term 
quadratic in q x and q y is set to zero since the curvature 
of the FS vanishes at the inflection point. The quartic 
term of q x controls the FS nesting. Hence, if we decrease 
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the coefficient b to zero, the CDW susceptibility will di- 
verge, and a phase transition to a CDW state will be 
accessed by tuning the parameter b. If we start from 
a nematic phase, b can be used to control the nematic- 



Here the constant 7 = (2/b) 1 ^/(2n). As we just men- 
tioned, 6 = implies a nested FS and this is the rea- 
son that 7 diverges as b — > 0. In contrast to the or- 
dinary Qs = 2kp incommensurate case studied before, 
there is no cancellation inside the particle-hole contin- 
uum for the inflection points, so that the naive scaling 
law, [q x ] = l,[q y ] = 3, and [u>] = 3, is always valid. 
The non-Gaussian terms now start at the order |$| 9 / 4 
as shown in Appendix [Cj This term is irrelevant at tree 
level. 

The low-energy heat capacity contributed by this mode 
has C ~ T 4 / 3 . The fermions near the inflection points 
have E" ~ |w| 13 / 12 , which means the fermionic quasipar- 
ticles near the inflection point will remain well defined, 
even when the nematic-smectic QCP is reached. 

VII. THE ELECTRONIC SMECTIC PHASE 

In the electronic smectic phase, i.e. in a conducting 
stripe phase, the electronic structure of the FL quasi- 
particles is changed by the modulation imposed by the 
smectic order parameter, and an energy gap will develop 
at the locus of the former FS. As a result new electronic 
bands along the direction of the ordering wave vector will 
be formed, with a gap in the electronic energy spectrum 
As oc gs&- As shown in Fig. [H^c), if 2kp/Qs is not 
close to an integer value, the highest band will have a 
closed FS, i.e. an electron (or hole) pocket, while the 
lower band (or bands) has instead an open FS (an "open 
orbit"). For 2kplQs close to an integer, as shown in 
Fig. Hfd), there will only be an open FS. The reconstruc- 
tion of the FS of the effective quasiparticlcs of a stripe 
phase in the cuprate high temp erature superconductors 
has been discussed recentlv |l43j as a possible explanation 
of the observation of Shubnikov-deHaas and deHaa s-van 
Alphen oscillations at relatively high magnetic fields jl44l |. 

We divide our discussion into two cases of interest: a) 
a smectic state in the continuum and b) a smectic state 



smectic transition. 

The Gaussian terms in the effective Lagrangian density 
of the smectic order parameter field is 



(6.15) 



(6.16) 



on a two-dimensional lattice. 



A. The electronic smectic phase in a continuum 

system 

We consider first a system with continuous rotational 
and translational symmetries which are partially and 
spontaneously broken by the smectic order parameter. 
We consider only the collective modes of the smectic or- 
der. The nematic collective modes will be added a poste- 
riori using symmetry considerations. In this section we 
derive the effective order-parameter theory for the elec- 
tronic smectic. The details of the calculation are pre- 
sented in Appendix iDl 

Let us consider a configuration of the complex smectic 
order-parameter field which is a small deformation of 
the ground-state configuration with expectation value "5, 

* = #[l + ^)e i **. (7.1) 

Here 5§ measures the amplitude fluctuations of $ and 
0$ measures the phase fluctuations. 

At the mean-field level, i.e. ignoring fluctuations, the 
main effect of the presence of a non-vanishing <E> is that 
the original FS becomes folded along the direction deter- 
mined by the ordering wave vector and a more complex 
electronic band structure results. The inter-band and the 
intra-band scatterings of the fermions will have different 
contributions to the effective low-energy theory. For the 
smectic Goldstone mode, 0$, in the static limit, the intra- 
band scatterings have no contributions but the inter- 
band scatterings will give a constant, which will cancel 
the constant coming from the terms As|$| 2 +its|<i>| 4 for 
Qs < 2k F , or A s |$| 2 + u s J$| 5 / 2 for Q s = 2k F , provided 
the expectation value of $ satisfies the self-consistency 
equation of the mean-field theory. This exact cancel- 
lation of the constant term for the Goldstone mode is 



which implies the smectic susceptibility 



%+ a -f+ b 4+^io ++ ^ qy+ f + b -f-.-m 
1 
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required by the Ward identity for translations along the 
ordering wave vector. The amplitude mode however will 
not have this cancellation, and a term proportional to 
5Q 2 will appear in its effective action. 

For u> ^ and q ^ 0, the intra-band and inter-band 

I 



where ip q is the angle between q and the stripe direc- 
tion and -/Vs(O) is the density of states in the smectic 
phase. A(ip q ), B(ip q ), and ns{<Pq) are coupling constants 
that depend on microscopic details and the direction of 
q, which reflects the anisotropic nature of the smectic 
phase. We will neglect the direction dependence of these 
three coefficients since they result in irrelevant contribu- 
tions at low energies and long distances. The first term 
oc u 2 in Eq. (|7.2[) is due to inter-band scattering, while 
the second, oc i\uj\q, is due to intra-band scattering. The 
intra-band contribution does not have the typical i\u)\/q 
form, because the contributions from the intra-band scat- 
tering vanish in the q = limit (see Appendix [Pi) The 
dynamic critical exponent is z = 2 at the QCP but z = 1 
in the smectic phase. This discontinuity is reflected in 
the singularity of Eq. (|7T2"]) at $ = 0. 

This behavior of the smectic Goldstone mode is very 
similar to that of the spin Goldstone mode in the 
commensurate antiferromagnetic phase (commensurate 
SDW) of a 2D FL studied by Sachdev, Chubukov, and 
Sokol [l45j] . In fact, although they are studying the 
commensurate SDW, while we are studying incommen- 
surate CDW, if the nematic fluctuations in our system 
are gapped by a lattice background or external fields 
which breaks the continuous rotational symmetry, the 
nematic-smectic transition we studied here shares many 
commo n fea tures with their work (especially "type B" 
in Ref. [145| |). Both these two transitions are driven by 
FL instabilities with a finite ordering wave vector and 
fermions provide similar non-analytic damping terms for 
low-energy bosonic excitations in both systems. For both 
cases, the QCP has z — 2 and the ordered phase has a 
Goldstone mode with z = 1. 

An RPA calculation (see also Appendix [D]) shows that 
the amplitude mode fluctuations of the electronic 
smectic order parameter has an effective low-energy La- 
grangian density of the form 

£ 5 * = g 2 s N s (0) (M - |A S |) \8${% co)\ 2 . (7.3) 

The dynamical term i\uj\/q comes from the intra-band 
scattering. In contrast with the case of the Goldstone 
mode 0$ , there is no cancellation of the intra-band con- 
tributions at q = for the amplitude mode 5$. Inter- 



scatterings also yield different contributions. The effec- 
tive low-energy Lagrangian density of the smectic Gold- 
stone mode (a real field in position space and time) 
takes the form (see Appendix |D|) 



(7.2) 

I 

band scattering contributes an irrelevant ~ lo 2 dynamical 
term which is ignored. The longitudinal CDW suscepti- 
bility in the smectic ordered phase is 

Xf (<f,w) = - i{5$(q,u)8$(-q,-u)} Tet 



9 2 S Ns(0) - |A S |) ' 

where q is measured from the ordering wave vector Qs 

Hence, in the electron smectic phase both the phase 
and amplitude modes of the smectic order parameter 
scale with a dynamical exponent z — 1. This result is 
interesting, because in most other examples of contin- 
uous symmetry breaking, the amplitude mode is either 
gapped or has a smaller dynamic critical exponent z than 
the Goldstone mode. As a result, at lower energies, the 
amplitude mode will become weaker, and can be ignored 
from the effective low-energy theory. However, in this 
case, since both the amplitude mode and the Goldstone 
mode have the same dynamic critical exponent, no mat- 
ter how low the energy scale is these two modes can no 
longer be separated. 

Nevertheless, in our problem, the amplitude mode is 
irrelevant for the following reasons. Using Eqs. (|5.5|) 
and (|7.1j) . we see that both the amplitude CDW fluctu- 
ations 5<E> and the CDW Goldstone modes 0$ couple to 
fermionic CDW operator n(q, u>) in the form: 

g s n{q, u) 6$(-q, ~uj) + h.c, 

ig s i> n(q,u>) <p<s,(-q, -to) + h.c, (7.5) 

but the coefficients in the action of <5$ (9sNs(0)) 
and_ 0$ {gs\$\Ns(0)A/k F , gs\$\N s (0)B /k F and 
gg\Q\ 2 Ns(0)ns) have different scaling behaviors. With 
z = 1, the coefficients for have dimension 3, but 
the coefficients 0$'s coefficients have dimension 1. As 
a result, the coupling between the smectic Goldstone 
mode 4><j> and the fermions is actually more relevant 
than the coupling to the amplitude mode 5$. Hence, we 
will not consider 5$ in what follows. 

We finally need to couple the smectic Goldstone mode, 
0$, with the nematic Goldstone mode 4>n- Similar to 
the classical case, 4>n couples to the electronic smectic 



k F gs 1*1 



iA(cp q ) 



k F gs\^\ 



Ks{Vq)<l \(j><s>{q,u)Y 
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field $ as a gauge field (Eq. (|4.15|) V Hence, in the elec- 
tron smectic phase, upon substituting Eq. (|7.1[) into Eq. 
(|4.15p . the nematic Goldstone mode 4>m develops a mass 
term ~ |$| 2 0^, i.e. it acquires a gap. As in the clas- 
sical case, the fluctuations of the nematic order produce 
bending of the smectic order, which costs energy. Thus, 
here too, this effect leads to the expulsion of the nematic 
Goldstone mode. 

Since the dynamics of the nematic Goldstone mode 



4>n has the structure i\co\/q, the mass term |$| 2 (£jy- of 
the nematic Goldstone mode in the smectic phase makes 
this mode to have z = 1 dynamics, just as the smectic 
Goldstone mode </>$ has. However, if we are close to the 
QCP, where the smectic amplitude fluctuation $ is small, 
the energy scale of the smectic Goldstone mode will 
be lower than that of 4>n . Hence, upon integrating out 
the nematic Goldstone mode 4>Ni we obtain an effective 
low-energy action for of the form 



g 2 s m 2 N s (0) B 



k F gs\<£ 



+ iA- 



KiK 2 q x + q^ 
k F g s \$>\ 
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(7.6) 



where A, B are two coefficients, and k\ and k,i are the 
defined in Eqs. (|4.17[) and (|4.18[) . Here, as in the classical 
case, the coefficient of the q 2 term vanishes by rotational 
invariance. 

By inspection of Eq. (|7.6p we see that the scaling di- 
mensions now are [q x ] — 1 and [q y ] — [oj] — 2. Hence, we 
can conclude that the fluctuations of 4><s> contributes to 
the low-temperature specific heat which scales with tem- 
perature as C ~ T 3 / 2 . It follows from this scaling analy- 
sis that in this case the system is above its upper critical 
dimension. Hence, true long-range order exists. The log- 
arithmic corrections g ener ated by the higher-order terms 
in classical smectics |l32j will not be present. Notice 
that these scaling dimensions are the same as what we 
got from the phenomenological theory in Sec. IIV[ before 
it was coupled to the fermions. However, the physics is 
now very different. Most important of all, in the pres- 
ence of fermions, the smectic Goldstone mode is damped, 
which is not the case in the phenomenological theory. 

A one-loop calculation of the fermion self-energy shows 
that the fcrmionic quasiparticles in the electronic smectic 
phase is a FL ( Appendix IE 2[) . Indeed, for most of the 
FS, the quasiparticle scattering rate at low frequency is 
E"(fcf,u;) ~ w 2 log|o;|. At the special points on the FS 
where the Fermi velocity is parallel to the ordering wave 
vector Qs {i.e. perpendicular to the stripes) ( marked in 
Fig. (He)), H"(k F ,uj) ~ M 3 / 2 for small u. However, this 
|w| 3 / 2 scaling will not be observed for 2}cf/Qs close to 
an integer, since no point on the FS will have vp || Qs 
as shown in Fig. [2jd). 



B. The electronic smectic phase on a lattice 

For most physical electronic systems the continu- 
ous rotational symmetry is reduced to a discrete point 
group symmetry of the underlying lattice. As it is well 
known [3, |33| the explicit breaking of rotational invari- 
ance by the lattice has important consequences for an 



electron nematic state. For the simple square lattice, the 
0(2)/Za symmetry of the continuum nematic state re- 
duces to an Ising-like Z2 symmetry, an "Ising nematic." 
On the other hand, at the nematic-FL QCP symme- 
try breaking by the underlying lattice are irrelevant (36| . 
provided the quantum phase transition remains contin- 
uous which, in many instances does not appear to be 
the case(37l 38]. Lattice effects manifest in the electronic 
structure and hence on the allowed shape of the FS of 
the FL state. Lattice anisotropics will then break the 
continuous rotational symmetry and act as explicit sym- 
metry breaking fields. In this case, the nematic phase 
does not have a true Goldstone boson which now be- 
comes gapped. At the level of the effective theory this 
effect shows up by the presence of a term proportional 
to <fi 2 N in the effective action of the nematic Goldstone 
mode, which is just an allowed mass term since there is 
no Ward identity to prevent it. This term makes 4>n & 
z = 1 mode, in contrary to z — 3 without lattice. Hence, 
the "pseudo- Goldstone" mode 4>n has now become irrel- 
evant at low energies. This changes a number of things. 
For one, the non-FL behavior of the nematic phase in the 
continuum is replaced by an anisotropic FL state with 
well-defined quasiparticles. Naturally if the effects of the 
lattice are weak enough, at sufficiently high energies (or 
temperatures) they can be neglected and the non-FL ef- 
fects become detectable above this crossover. On the 
other hand, as the quantum phase transition to the elec- 
tronic smectic phase is approached the gapped nematic 
Goldstone modes become irrelevant and the low-energy 
physics is dominated instead by the fluctuations of the 
smectic mode, which in the continuum case is a higher- 
energy excitation. In other terms, one now obtains a 
"fluctuating stripe" regime close enough to the nematic- 
smectic QCP. 

In the electronic smectic phase, the lattice anisotropy 
makes the nematic Goldstone mode irrelevant, which can 
now be neglected in the effective low-energy theory. As 
a result, we only need to consider the smectic Goldstone 
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mode to understand the low-energy physics of the smec- 
tic phase. Depending on whether the smectic (CDW) 
is pinned down by the lattice or not, the physics of the 
smectic Goldstone mode behaves very differently. In the 
absence of lattice pinning, the Lagrangian density shown 
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in Eq. (|7.2p will be the proper low-energy theory of the 
smectic Goldstone fluctuations, so the scaling behavior is 
fee] = fey] = M = 1 an d the transverse smectic suscep- 
tibility is 
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where q is measured from the ordering wave vector Qs- 
This low-energy behavior is different from the electronic 
smectic in the continuum. If the lattice effects are weak 
enough, we expect a crossover with increasing energy 
scales from the behavior of a smectic coupled to the lat- 
tice to one without this coupling. 

By following the same type of analysis used in the pre- 
vious subsections we can conclude that the heat capac- 
ity has a contribution due to the fluctuations 0$ with a 
T 2 temperature dependence. To the one-loop level, the 
fermionic quasiparticles at the FS will have a scattering 
rate a;) ~ w 2 log|w| for u> — > at most part of 

the FS, while £"(fcp ,u) ~ ui 2 for some special points on 
the FS discussed in Appendix IE 21 This observation is 
consistent with a FL behavior (Appendix IE 2[) . 

For a pinned CDW, the cancellation of the intra-band 
scattering for no longer holds. Therefore, the dynam- 
ical term is ~ *M/<7 with z = 1. The coupling between 
and the fermions is irrelevant in the pinned smectic 
phase, and the FL picture will hold. 

In general, there are two sources for pinning the CDW: 
lattice and impurities. We will consider the lattice pin- 
ning first. Lattice pinning is relevant for commensurate 
CDWs. For incommensurate CDWs, however, some lat- 
tice distortions will be required to pin the CDWs. A 
distorted lattice can pin down an incommensurate CDW 
only if the CDW order parameter is large enough. This 
is because the energy gain by pinning vanishes as |$| 
goes to 0. On the other hand, the energy needed to 
distort the lattice will not go to zero with decreasing 
|i>|. For most conventional CDW materials, the CDW 
ordering is very strong at low T so that incommensurate 
CDWs w i ll always be pinned down by distorted lattice 
[l05l li~4rl [l47l [l48j and unpinned CDW only appears at 
finite T above a phase transition [lOTj . 

However, for our systems, there is a second-order phase 
transition from nematic to smectic. When we are close to 
the nematic-smectic QCP, the smectic order parameter 
will be small, so that an unpinned incommensurate CDW 
phase is stable against a lattice background at least when 
it is close to this QCP. When we are far from the critical 
point, the smectic order parameter becomes large, so the 
energy gain by pinning may exceed the energy cost of 
distortion. Hence, a pinned CDW phase may form. 



Throughout this paper we ignore the effects of 
quenched disorder. Impurities, and more generally disor- 
der, affect strongly all electronic liquid crystal phases as 
they couple linearly to the their order parameters, lead- 
ing to the destruction of these ordered phases and to 
glassy- type states @. For the case of the CDW phases 
this problem was studied extensively in the 1980s, for 
which pinning is releva nt fo r d < 4 for systems with 
short-range interactions Il49l | l50j ] and d < 3 for long- 
range Coulomb interactions |l5l| . In the case of these 
quantum phase transitions the effects of quenched disor- 
der, even in the clean limit, are only poorly understood 
and we will not explore these problems in this work. 



VIII. FINITE-TEMPERATURE CROSSOVERS 
AND THERMAL PHASE TRANSITIONS 

We now discuss the effects of a finite-temperature 
on the electronic nematic and smectic phases and on 
their quantum phase transition. In Fig. [5J we present 
schematic phase diagrams for the three cases of interest: 
a) the isotropic case (no lattice), b) the system with a 
lattice background (a square lattice in this case) and an 
unpinned electron smectic (CDW) and in c) the case of 
a lattice with pinned CDW. 

In the nematic phase and at the nematic-smectic QCP, 
the nematic Goldstone modes have a finite temperature 
equal-imaginary time correlation function, i.e. the trans- 
verse nematic susceptibility, 



Xl{q) =N^(4> N (q,T)c} >N (-q,T)) let 



E 



0g 2 N N(O) ^ V£A S in 2 2 Vq + + K 2 q. 



(8.1) 

where we have set ks to 1 so [3 — 1/T and the sum 
runs over all bosonic Matsubara frequencies, uj n — litnT , 
where n £ Z. The susceptibility x±(q) is singular as q — > 
(due to the ui n = mode, where it takes the asymptotic 
form ~ T/q 2 . This behavior suggests that the nematic 
order parameter field has power-law correlations at finite 
T. This conclusion agrees with the classical theory of 
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FIG. 5: (online color) The schematic phase diagrams at finite 
T for (a) an isotropic system, (b) a system with a lattice back- 
ground, and (c) systems with a pinned CDW in the smectic 
phase. The horizontal axis is a control parameter that drives 
the quantum phase transitions: Fermi liquid — > nematic — > 
smectic at T — 0. The vertical axis is the temperature T. The 
thick (blue) lines are phase transitions belong to the Ising (or 
q-state Potts model) universality class while thin (red) lines 
are KT transition phase boundaries. The dashed lines mark 
the crossover to the quantum critical regime (QC). In panel 
(a) we show that the nematic phase (in the absence of lat- 
tice symmetry breaking) only has long range order at T = 0, 
and it is a critical KT phase at all temperatures below the 
KT transition. In panel b) the lattice reduces the symmetry 
to Z2 and there is long-range Ising nematic order at finite 
temperatures. In this case there is also a finite-temperature 
smectic phase which is critical (KT) if the smectic (CDW) 
order is unpinned and has long-range order in the latter case 
(shown in panel c). 



a 2D nematic phase, which belon gs to the Kosterlitz- 
Thouless (KT) universality class (El ll52T |. 

By comparing the T/q 2 term with the equal-time cor- 
relation function at T — 0, which is ~ qlog(Ep /<?) where 
Ep is the Fermi energy, which enters the calculation as a 
high energy cutoff, we see that T = behavior becomes 
dominant when q 3> T 1 / 3 . 

From the theory of classical liquid crystals [J, [||, we 
know that in a fully translationally invariant system, 
e.g. in the absence of a lattice, there is no stable finite- 
temperature smectic phase in two space dimensions jl52| . 
which implies that the equal-time correlation function of 
the smectic order parameter decays exponentially fast as 
a function of distance at any finite T. Due to the effects 
of the proliferating dislocations of the smectic the actual 
finite-temperature phase is a nematic. 

In the nematic phase and at the nematic-smectic QCP, 
the finite T equal-time correlation function of the smectic 
fields, i.e. the smectic susceptibility, is 

X S (q)^(q,TM-q,T)) Ict 

= /3 ^ Co M + C x q\ + C y q 2 + A (8 - 2) 

Here, q is measured from the ordering wave vector Qs 
and A = A5 + /(T), where f{T) is a function of tem- 
perature which vanishes at T = 0. The leading term 
in f(T) comes from the quartic term of the smectic field 
~ usT and the gauge-like couplings between the nematic 
Goldstone mode and the smectic fields ~ Q^T, where the 
ordering wave vector Qs represents the strength of the 
gauge-like coupling. The absence of a finite-T smectic- 
ordered phase implies that A 7^ for T > 0. Hence, 
on the nematic side, at low enough temperatures A > 
and it is essentially equal to Ag > 0. However, in the 
quantum critical region (denoted by QC in Fig. EJ), the 
singular behavior of x S {%) is regulated purely by finite 
temperature. Notice that /(T) depends on irrelevant op- 
erators of the QCP; we conclude that the infrared di- 
vergence is regulated by irrelevant operators, including 
the gauge-like coupling between smectic and the nematic 
Goldstone mode, which implies that this gauge-like cou- 
pling is a dangerous irrelevant term near the QCP and be- 
come relevant at finite T. This conclusion agrees with the 
theory of the classical nematic-smectic transition at finite 
T, where the gauge coupling is known to be relevant. The 
boundary between the nematic and the quantum critical 
regime is determined by the condition f(T) ~ A5, which 
isT~ A s . 

On the smectic side, the equal-time correlation func- 
tion of the smectic Goldstone mode fluctuations, i.e. the 
transverse smectic susceptibility X±(o) (where q is the 
momentum measured from the ordering wave vector Qs 
of the smectic phase) is 
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X b ± (q) = m 2 (<Ps(q,T)<P s (-q,T)} Ict = 



A- 



(8.3) 



Kiq x 



The most divergent term in x±(o) has the asymptotic 
behavior T j {niq^ + ^q 2 ), which implies that the Fourier 
transform of x±(<l)> the transverse susceptibility of the 
smectic Goldstone mode, is infrared divergent at finite 
T. This divergence implies that the equal-time cor- 
relation functions of the smectic order parameter field 
($(x)t$(y)) exactly vanishes, for all f / i/. Thus, this 
system exhibits a form of "local quantum criticality" , i.e. 
where the correlation length of the equal-time correla- 
tion function essentially vanishes but the equal-position 
auto-correlation function scales as a function of time. In 
Fourier space, this means that the finite temperature dy- 
namical susceptibility scales in frequency but not in mo- 
mentum. This behavior was discussed recently in the 
context of the FL-nematic QCP[!i| and in the quantum 
critical behavior of quantum dimer models jl53l |. 

The vanishing of the equal-time correlation length 
however invalidates the assumption of small smectic 
Goldstone fluctuations, where topological defects (dis- 
locations) become important. With dislocations, the 
smectic order parameter vanishes (as translation invari- 
ance has been restored) and the equal-time correla- 
tion functions of the smectic order parameter become 
short ranged. This physics proceeds, as usual, by the 
non-perturbative Kosterlitz-Thouless mechanism of de- 
fect proliferation. At low temperatures, where the den- 
sity of the dislocations is low, the correlation length di- 
verges, which recovers the proper behaviors of the zero- 
temperature quantum smectic phase. 

For systems with a lattice background, the nature of 
the thermal nematic phase transition is determined by 
the point group symmetry of the lattice, e.g. for the 
square lattice it is an Ising transition. Also in the pres- 
ence of a lattice there is a stable electronic smectic phase 
at finite temperature. In this case the finite-temperature 
nematic-smectic transition is in the KT universality class 
for an unpinned CDW. With a lattice background, the 2D 
unpinned smectic will be at its lower critical dimension. 
Therefore, the smectic phase has quasi long-range order 
and a finite-temperature KT transition is expected. For 
smectics pinned by lattice, the transition instead belongs 
to the universality class of the g-state Potts model, with q 
equa ls to n for a lattice with n-fold rotational symmetry 
|lldj |. On the nematic side, Eq. (|8.2[) is still valid, but 
A — > as we approach the finite-temperature phase tran- 
sition. On the unpinned smectic side, the most singular 
term of the correlation function of its Goldstone mode 
will be ~ T/q 2 , which now will be unaffected by topolog- 
ical defects in the finite-T smectic phase since in the case 
of a lattice the dislocations are confined. This behav- 



ior of a phase smectic phase mode leads, as usual, to an 
equal-time correlation functions of the order parameter 
with power-law behavior as a function of distance. In the 
high-temperature phase, where the dislocations become 
deconfined, the order parameter has only short-ranged 
spatial correlations. 

IX. DISCUSSION AND CONCLUSIONS 

In this paper we presented a phenomenological theory 
of quantum phase transitions in electronic liquid crys- 
tal phases. For simplicity we have only considered the 
charge channel and discussed the behavior of two charge- 
ordered phases: the electron nematic and the electron 
smectic. The latter phase has the same symmetries as a 
charge stripe phase and a unidirectional CDW. We have 
not discussed the behavior of spin excitations and the as- 
sociated magnetic phases, i.e. a spin stripe. The main 
results are summarized in Table HI 

We discussed in detail how to describe a fluctuating 
charge stripe phase, a metallic electron nematic phase 
close to a quantum phase transition to an electronic smec- 
tic (or stripe) phase. We analyzed the nature of the fluc- 
tuations of the bosonic collective modes in the nematic 
as this quantum phase transition is approached: it is de- 
scribed by a low-energy Goldstone nematic mode (with 
dynamic critical exponent z — 3) and a higher energy 
smectic collective mode (with dynamic critical exponent 
z = 2), and discussed in detail the behavior of the ne- 
matic and smectic correlation functions in this regime. 
We also gave a detailed analysis of the electron smectic 
phase and of the behavior of the correlation functions 
in this phase. In particular we give explicit form of the 
dynamical susceptibilities at zero temperature, includ- 
ing the dynamics induced by the fermionic fluctuations. 
This analysis was also done in the electron smectic phase. 
The behavior of the nematic and smectic susceptibilities 
should be useful to interpret experiments that can probe 
this dynamics, particularly light-scattering experiments. 
With some minor changes the form of these correlation 
functions and susceptibilities also apply to the analysis 
of magnetic fluctuations, such as neutron scattering ex- 
periments in the fluctuating stripe regime. 

We also developed a description of the quantum critical 
behavior at this phase transition to an electronic smec- 
tic phase. This effective critical theory is a quantum- 
mechanical generalization of the classical McMillan- 
deGennes theory for the nematic-smectic liquid crystals 
to quantum-mechanical metallic phases with the same 



25 



pattern of symmetry breaking. It is a theory at zero 
temperature and it includes the effects of strong quan- 
tum fluctuations which turn out to have a very different 
character than their classical counterpart. Due to the 
non-trivial effects of the fermionic fluctuations, this quan- 
tum phase transition is continuous whereas the classical 
finite-temperature transition is weakly first order. 

We presented an extensive discussion of the fate of the 
fermionic quasiparticlcs in each phase and at the quan- 
tum phase transition, and of the resulting non-FL behav- 
iors. The resulting non-FL effects are quite rich. These 
results were obtained within a perturbative expansion 
in powers of the coupling between the fermions and the 
order-parameter fluctuations. As such, the obtained non- 
FL behaviors represent primarily a breakdown of pertur- 
bation theory rather than asymptotically exact results. 
The reliability of these low-order results will be checked 
in the near future using non-perturbative approaches 
such as higher dimensional bosonization. In any case 
our perturbative analysis of the quasiparticle self-energy 
shows once again that electronic liquid crystal phases are 
naturally compatible with non-FL behavior. 

Although much of the theory that we presented is de- 
veloped in the context of a continuum system, i.e. a sys- 
tem in which the effects of the coupling of the underlying 
lattice on the electronic order is ignored, we also gave 
a qualitative analysis of the symmetry-breaking effects 
resulting from the coupling to the lattice. This analysis 
will be generally correct even though the details of the 
microscopic band structures were ignored, as they will be 
reflected in the form of symmetry breaking fields and in 
their coupling constants. We expect that the discussion 
of the continuum (isotropic) theory will be applicable if 
the effects of the coupling to the lattice are comparatively 
weak (as expected far from van Hove singularities) and 
at temperatures high compared to the scale of these cou- 
plings. Quite surprisingly, although the pinning of the 
nematic Goldstone mode by the lattice was expected, we 
found that the smectic collective modes in general re- 
main strongly fluctuating as the QCP is approached and 
into the smectic phase. This happens provided that the 
magnitude of the ordering wave vector obeys Qs < 2k,F- 
Contrary to what happens in the case of the FL-nematic 
quantum phase transition, in which lattice effects often 
drive the transition first order, the nematic-smectic tran- 
sition can naturally be continuous even in the presence 
of the coupling to the lattice. 

We also considered the effects of low-temperature ther- 
mal fluctuations, and showed that the signatures of the 
electronic smectic phase can be detected through heat- 
capacity measurements as subleading corrections to FL. 
In the electronic nematic phase and at the nematic- 
smectic QCP, the heat capacity is dominated by ne- 
matic fluctuations, which goes as T 2 / 3 . At the QCP, 
the smectic mode yields an additional contribution ~ T 
for Qs < 2kF, a subleading term ~ J 13 / 2 for commen- 
surate CDW with Q s = 2k F and ~ T 4 / 3 for the case 
of inflection points. In the electronic smectic phase, the 



smectic Goldstone mode gives a contribution to the heat 
capacity of C ~ T 3 / 2 in the absence of lattice symmetry 
breaking effects and T 2 with lattice symmetry breaking 
effects. 

The approach that we followed is semi- 
phenomenological and it is based on the picture of 
a Fermi liquid that is coupled to an effective field theory 
describing the fluctuations of the nematic and smectic 
(stripe) order parameters. This line of attack can be 
justified, at the level of mean-field theory, in weakly 
coupled systems based on the existence of a Fermi 
liquid for some range of parameters and its subsequent 
instabilities. Examples of this approach is the work 
on the FL/electron-nematic quantum phase transition 
of Oganesyan et al [33| in continuum models and of 
Metzner and coworkers [34 Hf| HI) HI] and H.-Y. Kee 
and coworkers (37i l38j in lattice models. The extension 
of these works to the electron nematic/stripe phase 
transition that we discuss here is possible and we have 
obtained some unpublished results in this direction. 
However, this quantum phase transition requires that 
a coupling constant be larger than a critical value, 
which typically is not small, and hence the reliability of 
mean-field methods in this regime is at least problematic. 
This will be discussed in a separate publication. 

We have not discussed the realization of these elec- 
tronic liquid crystal phases, and of the quantum 
phase transitions we discussed in microscopic models 
of strongly correlated systems. Nevertheless sufficiently 
close to a continuous quantum phase transition this ap- 
proach is likely to give the correct universal behavior. 
However, the use of theories based on the breakdown 
of the FL state is problematic in the strong correlation 
regime. In addition, so far there is no microscopic model 
in which both nematic and stripe phases are known to 
occur. The existence of an electron nematic phase in a 
microscopic model of a strongly correlated system has 
only been shown for the strong coupling regime of the 
Emery model of the cuprates in the asymptotically low 
doping regime [4J|. The existence of a metallic stripe 
phase in the same model for the doping range x ~ 1/8 
has been suggested by a variational wave function (and 
hence mean-field in spirit) approac h wh ich projects out 
double occupancies of Cu sites 0, Il54j . This suggests 
that the strongly coupled Emery model may more gener- 
ally exhibit both nematic and stripe phases in its phase 
diagram. A number of publications have attempted to 
describe both stripe and nematic phases in Hubbard, 
t — J, and Em ery models using slave-particle methods 
[I55[|l56lll57l ]. However, in spite of their widespread use 
in the field, slave-particle mean-field theories are notori- 
ously unreliable. 

It would be highly desirable to have high quality nu- 
merical simulations to address this problem in models of 
strongly correlated systems. Density matrix renormal- 
ization group (DMRG) calculations have provided strong 
evidence for stripe correlations in Hubb ard-t ype models 
on narrow strips (with up to 5 legs) (l3l. Il58| . However, 
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the geometry used in DMRG, which breaks the rotational 
invariance under 90° rotations of the square lattice ex- 
plicitly makes it difficult to distinguish a stripe from a 
nematic phase. The same problem arises in finite-size di- 
agonalizations of small systems. Quantum Monte Carlo 
simulations are less affected by such geometric limitations 
but suffer from the notorious fermion sign problem at low 
temperatures. QMC simulations have indeed shown an 
increase in nematic fluctuations at low temperatures in 
Hubbard type models (see a discussion in Ref. [l3j]) but, 
as far as we know, not yet in the Emery model. 
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APPENDIX A: TENSOR FORM OF THE ORDER 
PARAMETERS 

In this appendix we rewrite the order parameters in a 
tensorial form which makes their correct symmetry trans- 
formation properties apparent. Much of what is done 

I 

f dfcdfi f dqdtJ 1 / / k 2 x k x k y 



Again, it is easy to check that only the traceless part of 
the tensor composed by the momenta k is needed here, 
since the trace of that tensor has no contribution. Notice 
that although there is a factor of 1/k 2 in our formula, it 
will not cause divergence, since what we are interested 
in is a CDW, a state that orders at a finite wave vector 
k~ Q s . 

Finally, using the fermionic density quadrupole tensor 
defined in Eq. (|5.ip , Eq. (|5.3[) can be written in tensorial 
form as 



-g N / drdt tr[QN] (A4) 



here follows closely the analysis of the classical case[5[. 

In 2D, because the rotation group 50(2) group is 
Abelian, and has only two ID representations with 1=2 
(l z = ±2), it follows that the nematic order parameter 
N, a 2 x 2 symmetric traceless tensor, has just two in- 
dependent components, n\x and ni2, as shown in Eq. 
(|4.1|) . The Abelian nature of the 50(2) group enables us 
to use complex numbers, instead of tensors to represent 
the action of the group, as shown in Eq. (|4.2p . 

All the formulas we have presented in the main text us- 
ing the complex order parameter can be translated into 
the tensor language, which can more naturally be gener- 
alized to higher dimensions. Thus, Eq. (|4.7[) becomes 

-2k tr {N(f) [DN(f)N(r)]} . (Al) 

D is the rank 2 tensor formed by the outer product of 
the two-component real vector (d x , d y ): 

D=(^)®(0*A), (A2) 

In Eq (|A1|) the first derivative operator acts on the second 
N factor, while the second derivative operator acts on 
the last N factor. Here only the traceless part of D, 
D — tr[D]I/2, gives non-zero contribution to Eq. (|A1[) . 
The assignment of derivatives is not unique. But only 
this term gives a linearly independent contribution to the 
nematic Goldstone mode. Other assignments either have 
no contribution, or just give a contribution proportional 
to this one. Therefore, there is no needs to study other 
terms. 

Similarly, Eq. (|4.12[) has the tensor form: 



H p(k - q/2, Q - u/2)p(-k- q/2, -SI - ui/2). (A3) 
I 

APPENDIX B: QUANTUM CRITICAL POINT 

FOR Q s < 2k F 

As discussed in the main text, at the nematic-smectic 
critical point Qs < 2kp, the nematic Goldstone mode is 
the low-energy mode. To obtain its effective low-energy 
effective action we integrate the high energy mode of the 
smectic field, $. Terms like 

{{d x & d x <5>))^ N( j) N , (bi) 

will generate the leading corrections to the <pN propaga- 
tor. We evaluate the integral numerically to deduce the 
kernel for the quadratic fluctuations. To the one-loop 
level, it is 

C iq 2 + iC 2 ^\\og{u J /q 2 )\/q\ (B2) 
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for < oj <C q 2 ■ The first term C\q 2 renormalizes the 
two Frank constant K\ and K 2 , while the second one is 
subleading compared to the dynamical term of the ne- 
matic Goldstone mode, which is proportional to i|w|/g. 
This result suggests that the smectic field $ is an irrele- 
vant perturbation to the nematic Goldstone mode at the 
nematic-smectic critical point. 
Terms of the form 

(ft,**^**^)*** (B3) 

provide leading order corrections to the propagator of the 
smectic field $, and yield a self-energy correction to the 
field <f>. The loop integral is computed numerically and 
is well fit by the form 

~C 3 q 2 + iCiLU 5 / 4 , (B4) 



when w>5 3 . The first term renormalizes the constants 
C x and C y defined in Eq. (|4. 15|) and the second one is 
subleading compared to the dynamical term of the smec- 
tic field which is linear in |w|. 

Therefore, these arguments provide strong evidence 
that the coupling between the nematic Goldstone mode 
</>jv and the smectic order parameter $ is irrelevant. In 
the low-energy theory, we can treat in practice (f>N and 
<& as two separate modes. 



APPENDIX C: NON-ANALYTIC TERMS OF 
THE EFFECTIVE FIELD THEORY OF THE 
ELECTRON SMECTIC WITH Q s = 2k F 



The mean-field Hamiltonian of the smectic phase is 

H MF = J-^{e{k)^(k^(k)+g s § ^(k+QsMk) + gs ^(fc-tfsM*)) + A s |#| 2 . (CI) 

The smectic order parameter $ will act as a periodic background potential in the direction of Qs which will reconstruct 
the band structure and the FS. For Qs — 2kp, we will only consider the lowest two bands. By ignoring higher bands, 
we can diagonalize this mean-field Hamiltonian to get the single particle dispersion relation 



E±(k) = - (e(k) + e(k + Q S )) ±^j- (e(k) - e(k + Q s )) + 9 2 s m 2 - (C2) 



Here the + sign is for the upper band and the — sign 
is for the lower band. The lower band is partially filled 
by the fermions, but the upper band is empty. So the 
Landau free energy will be 



F(\<b\) = 



E_(/c)<M 



d 2 k 
(27T) 2 



£_(fc) + A s |$| 



(C3) 



where /i is the chemical potential. By expanding the 
dispersion relation around the two points marked on Fig. 
Mb) 



eiAo) = ±( ly + Till, 



(C4) 



where q is measured from the special points marked in 
Fig. [Hb). The Fermi velocities at these two points are 
just opposite to each other, vp = ±1, and the curvatures 
of the FS have the same value k at these two points. 

By using this approximate dispersion relation, the Lan- 
dau free energy can be determined as 



F(|$|) = F(0)+r|*| 



|$|6/a. 



(C5) 



Here the coefficient of the quadratic term r depends on 
the high-energy cutoff, which reflects that the Landau 
free energy depends on the band structure all the way 
down to the bottom of the band. 



On the other hand, v! is universal, as it depends only 
on the dispersion relation around the two special points, 



r(-i/4) 



57r 3 / 2 r(i/4) 




0.049 



5/2 

9s 



(C6) 



u' diverges in the limit of k — > 0, which means terms in 
lower order, such as |$| 9 /4, will be generated (see below). 

Around an inflection point, the dispersion relation is 
shown in Eq. (|6.f4p . Similar to the calculation above, a 
term proportional to |$| 9 / 4 will be found. The coefficient 
of these term u" will be 

2r(- 



-1/8) 9 T 



97r 3 / 2 r(3/8) 6V4 



9/4 

0.15^ (C7) 



APPENDIX D: GOLDSTONE MODE IN THE 
SMECTIC PHASE 

In the smectic phase, the smectic order parameter ob- 
tains an expectation value. We write $ as 



$ = $ + <5$. 



(Dl) 



$ can be considered as a periodic potential background 
for the fermions. As a result, the fermions will form band 
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structure in the y direction. Define the Bloch states as 
Mk) = J2 T n , m (k)e m U(k + mQ s ), (D2) 

m 

where T n . m is an orthogonal transfer matrix which de- 
pends on the amplitude of the order parameter, and <f> is 
the phase of $. In general, the Bloch wave is related to 
the plane wave by a unitary transformation. Here, due to 
the fact that only one harmonic of the CDW with wave 
vector Qs is considered, after a proper spatial transla- 
tion, which is a shift of <j>, it can be simplified to an or- 



thogonal transformation, which enables us to define the 
orthogonal transfer matrix T n m . The inverse formula 
can also be written down, since the T n ^ m e" n ^ must be 
unitary: 

i>{k + mQ s ) =^T n , m (fe)e- iro ^n(fc)- (D3) 



Integrating out the fermionic degrees of freedom, the 
leading term in the power series of the smectic fluctua- 
tions <5$ starts from the quadratic order. 



V ^ ' Kq U \-n'(-Q s -q,-oj) A s -U(Q S -q,-uj) J \8&(-q,-u>) J K 1 

Notice that 5&(q,u) ^ 5$(-q,-w), since the field $ is complex. The continuous translational symmetries have 
been broken into discrete ones. Therefore, a process which changes the momentum by nQs, where n is an integer, is 
allowed. This is the reason why we have the terms such as 6<&'(q ) cj)6$*(—q, —u>). 

n(Qs + <f; in the diagonal terms is the standard fermion bubble integral in the smectic phase but with momentum 
close to Qs, 

[dkxdCLt [dk 2 dU 2 

(V (fci + Qs + q/2, Qi + w/2)i>{k x - q/2, f*i - u>/2)^(k 2 - q/2, il 2 - u/2)^{k 2 + Q S + q/2, il 2 + lj/2)) . 

(D5) 

Since the eigenstates are Bloch waves, we need to transfer tjj into ijj n . 
U{Q s + q.u) 

/ Jj^ T ni, mi +i{k + q/2)T n2 ^ mi {k-q/2)T n2 . m2 {k (D6) 



mi ,m2 ,ni ,n 2 



where F ni ^ n2 {k,q,uj) describe the scattering between fermions in band n\ and n 2 , which is defined as 



2tt 



1>U (* + <f/2, n + w/2) ip ni (k + q/2, ft + w/2) ) (v>I 2 (k - q/2, Q - w/2) ^„ 2 (k - q/2, Q - uj/2 



(D7) 

The off-diagonal H'(Qs + q, oj) is very similar to Ii(Qs + <f> 

^ /-dfexdOi fdk 2 dn 2 

n(QsHu)= iwiw 

<V (*i + Qs + 9/2, fii + w/2)V(fci - 9/2, fii - uu/2)^(k 2 +Q S ~ q/2, il 2 - w/2)^{k 2 + q/2, ft 2 + lu/2)) . 

(D8) 

Again, we need to transfer ip into ip n . 

n'(<5s + f,w) 



/* d/c 

X] / 7^32 T ni,mi+l(£ + 9/2)T n2 ^ rni (k - (f/2)T„ 2 , m2+ i(fc - q/2)T num2 (k + q/2)F ni , n2 {k,q,w) 



mi ,7712 )7l2 



(D9) 



To study the low-energy excitations, we need to diago- 



nalize the matrix of Eq. (|D4[) . Then two eigenmodes ap- 
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pear. In the limit where we take ui — > first, then q — > 0, 
it is easy to see that eigenvalues of these two eigenmodes 
are U(Q S , 0) + A s ± \W(Q S , 0)|. The equation 



n(Q s ,0) + As-n'(Q s ,0) = 



(D10) 



reproduces the mean-field self-consistent equation which 
determines <&. Therefore, the mode which takes the mi- 
nus sign will be the Goldstone mode and the other is the 



J 



amplitude mode. 

As for the Goldstone mode, the intra-band scattering 
-F ra ,n and the inter-band scattering F nitTl2 where n\ ^ ?i2, 
have different contributions. The intra-band scattering 
have no contribution in the q — ► limit, because its con- 
tribution to H(Qs,0) and |n'(Qs,0)| cancels. 

When n\ = ri2, term in H(Qs, 0) is 



dk 



^ ] (27r) 2 ^™> mi +i (^"i"'! (^)-^t»,fn2 {k)T n rn2 -i r i(k)F n (fc, 0, 0) , 



(Dll) 



mi ,7712 ,Ti 



and term in II' (Qs, 0) is 



E 



dfe 



(k)T nymi (k)T n:in2+ i (k)T n , m2 (k)F n (k, 0, 0). 

I 



(D12) 



From these two equations, the contributions are canceled 
for the Goldstone mode. 

The inter-band scatterings are gapped, so that we have 

n f (E ni (k))-n f (E n2 (k)) 

F„ llBa (M,0) = V - 1 - >, S (D13) 

£ ni (fe) - £„ 2 (fc) 

where E n (k) is the eigen-energy of fermions in band 
n with momentum k and n/ is the Fermi distribution 
function. Summing all the inter-band contributions, the 
equation H(Qs, 0) + As — H'(Qs, 0) = reproduces the 
mean-field self-consistent equation which determines 5> 



J 



£ = g 2 \<b\ 2 N s (0) [B(<p 



as required by the Ward identity 

Now let us study the frequency dependence of the 
smectic Goldstone mode. It is straight forward to show 
that the frequency dependence from the inter-band scat- 
terings will start from the order oj 2 because the inter- 
band scattering has an energy gap. 

We expect the intra-band scattering gives a dynamical 
term ~ i\ui\/q to II and II' since there is no energy gap 
and the T matrix has no singular points. However, as we 
just showed, the intra-band scattering contributions to II 
and II' cancels at q = 0. Hence, the leading dynamical 
term is ~ i|cj|g. The Lagrangian density reads 



k F gs |$| 



iA((p q 



\LU\q 



k F gs |$| 



(D14) 



Here $ is the expectation value of the smectic order 
parameter and the smectic Goldstone mode <f>^ is defined 
in Eq. (|7.ip . (p q is the angle between q and the stripe 
direction and Ns(0) is the density of states in the smectic 
phase. A(ip q ), B(ip q ) and Ks(<p q ) are coupling constants 
that depend on microscopic details and the direction of q 
which reflects the anisotropic nature of the smectic phase. 

This low-energy theory has z = 1 which is different 
from z = 2 at the QCP. Notice that the denominator 
of the first two terms in the above formula contains |$|, 
which just means that the z = 1 breaks down when on 
approach the QCP from the smectic side. On the other 



side, if we notice that there is a coefficient |<I>| 2 in the front 
of the Lagrangian density, the |$| in the denominator 
cause no divergence at small |$|. In fact, as we can see 
that the intra-band scattering vanishes as |$| goes to 0, 
which is what one should expected. 

It is easy to check that the coefficients of A(tp q ) and 
B(ip q ) are dimensionless, while B(tp q ) and Ks(<fq) are 
real. A(tp q ) is a real function for most values of ip q . 
However, from some special cases, for example Fig. [2jd) , 
there will be no FS whose Fermi velocity is close to the y 
direction. As a result, the particle- hole excitations with 
momentum close to the x direction will not be damped. 
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Therefore, A(<p q ) will be imaginary for (p q ~ or 7T. 



APPENDIX E: RPA CALCULATION OF THE 
FERMION SELF-ENERGY 



scribed by the Landau FL theory. But if n is less than 1, 
the self-energy correction will dominate and the Landau 
FL theory will break down at low energies. As a result, 
the system will exhibit non-FL behavior. 



The imaginary part of the fermion self-energy correc- 
tions at the FS, H"{kp,us), scales with frequency (as 
us — > 0) as ~ M^. When the scaling index [i is larger 
than 1, the low-energy theory of the fermions can be de- 



1. The nematic-smectic QCP 

To the one-loop level, for us > 0,the imaginary part of 
the fermion self-energy correction from $ is 



J 



E s (fc F ,c) = ^ 



0<e(k F +Q s -q)<ui 



dq 



r 



B s (q,u> - elk 



(El) 



where Bs(q, us) is the spectral density of the smectic field. 
At the nematic-smectic QCP (Qs < 2kp) Eq. (|6.8[) gives 



Bs (q, w) 



2C Q us 



(G x ql + C y qlf 



(E2) 



we find that for the 
Eg - \us\ 1 / 2 , while for 
2 For the special case 



By substituting this to Eq. (|ET 
four points marked on Fig. [2Ja 
all other points on the FS So ~ us 
of C x = C y — C, an analytical formula can be achieved. 
For the four special points, Eg = y / 2|Z7f/ (Att^/CCq) valid 
for \us\C/C < 1. For other points on the FS, E'^jfep) = 
us 2 v 2 Cq I \1&ttG 2 A 3 ) , where the constants A and v comes 
from the expansion of the dispersion relation near kp + 
Qs as e(q + k F + Q s ) = M + A + • q + . . .. 

The nematic Goldstone mode has an us 2 / 3, contribution 
to E" for most points on the FS except the four points 
in the main-axis directions where E" 



M 2 / 3 behavior is dominant over a Eg 



us\ 3 / 2 [33]. The 
us 2 scaling, but 



not over a Eg 



us] 1 / 2 contribution. As a result, the four 
main-axis directions still have FL behavior. The special 
points shown in Fig. [2fa) will have E" ~ \us\ x / 2 non-FL 
behaviors. For all other part of the Fermi surface, E" ~ 
M 2 / 3 due to the nematic Goldstone mode. If the nematic 
Goldstone mode is gapped due to a lattice background or 
an external field, the E" ~ |u;| 2 / 3 non-FL behavior will 
disappear. But the E" ~ l^j 1 / 2 non-FL behavior at the 



four special points will persist. 

The incommensurate CDW of Qs — 2kp case has a 
relevant perturbation <£> 5 / 2 , so the critical theory will be 
controlled by microscopic details. RPA calculation will 
not be reliable in such a situation. 

However we know that the commensurate CDW with 
Qs = 2kp is a second order transition described by the 
Gaussian theory, so RPA is applicable. Here we follow 
the same approach mentioned above. For most part of 
the FS, this bosonic mode will provide a self-energy cor- 



rection Ec 



for small us. However for the special 



Eg is linear in us when usn -C 1- \us 



2/3 



point k = ±kpe y , 
term from the 4>n mode will not be present for the com- 
mensurate CDW due to the lattice background. Hence, 
the fermions at this critical point are described by a FL 
theory, except for fermions with vp \\ Qs- 

For the special case of the inflection points, by follow- 
ing the same procedures, we get E" ~ |w| 13 / 12 at inflec- 
tion points, where the FL picture is valid for the whole 
FS. 



2. The electron smectic phase 

In the smectic phase, for us > 0, the fluctuations of 
the smectic Goldstone boson will contribute to the 
fermion self-energy: 



J 



1*1 



0<<L n i(k-q)<ui 



dq 

(27T) 2 



9s' n ( fc ;<?) 



B d 



q, \u\ - e n i{k - q) 



(E3) 



Here B^^us) is the spectral density function of the tf>$ mode presented below, e n (k) is the dispersion relation for 
fermions in band n and the vertex is 

5g™ (k,q) =ig s y^ y \r n ^ m+ i(k - q)T n ^ m (k) - T n>m {k - q)T n < , m+ i(k) 
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with T ni7n being the matrix element of the orthogonal transformation defined in Eq. (|D3|) 

The energy gap between two different energy bands dictated that the contributions from the interband scatterings 
to E" (n ^ n') scale as uJ 2 at small u>, which is subleading to the FL behavior. As for the intraband scatterings, the 
vertex vanishes linearly in the long wavelength limit (q ~ 0) as 



9s 



\k,q) =iJ2 [vT„, m (fc)T„, m+1 (fc) - VT n , m+1 (fc)T„, m (fc)] ■ q + 0(q 2 ). (E4) 



This is because the Goldstone mode cannot couple directly to the fermion density but to its fluctuations, as required 
by the translational symmetry. As a consequence of this structure of the vertex, we found that the quasiparticle 
scattering rate from intraband scatterings are also subleading corrections to the FL behavior by numerically evaluating 
the integral in Eq. (|E3p . 

For an electron smectic without a lattice background, Eq. (|7.6p gives 



B k F 2Aw y / Kl K 2 1 q*+ql 



For most part of the FS, the intraband scatterings lead to a fermion scattering rate E^ ~ w 2 log|w| in the small 
frequency limit. However for some special points, where the Fermi velocity is perpendicular to the stripe direction, 
E^ # ~ |w| 3 / 2 . This result suggests that the fermions will be described by a FL. For unpinned smectics with a lattice 
background, we have 

B (ff u )= kp 2A(<P q )ujq 

* ' 9smN s (0) (A^ q )toq) 2 +(B(<j) q )ui-k F gsmK S ((b q )q 2 ) 2 ' 



and that lim w ^oE"(a>) ~ w 2 log|ci;| at most part of the 
Fermi surface, while E"(/cf, ui) ~ uj 2 at the special points 
where g^' n (k,q) is independent of qt to the linear order 



of q, with q t being the component of q parallel to the 
reconstructed Fermi surface. 
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